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Abstract 

We investigate the possibility of electron pairing in small metallic nanoparticles at zero tem- 
perature. In these particles both electrons and phonons are mesoscopic, i.e. modified by the 
nanoparticle's finite size. The electrons, the phonons, and their interaction are described within 
the framework of a simplified model. The effective electron-electron interaction is derived from 
the underlying electron-phonon interaction. The effect of both effective interaction and Coulomb 
interaction on the electronic spectrum is evaluated. Results are presented for aluminum, zinc 
and potassium nanoparticles containing a few hundred atoms. We find that a large portion of 
the aluminum and zinc particles exhibit modifications in their electronic spectrum due to pairing 
correlations, while pairing correlations are not present in the potassium particles. 

PACS numbers: 
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I. INTRODUCTION 



Pairing effects in metallic nanoparticles have been investigated extensively since the ex- 
perimental work by Black, Ralph, and Tinkharoii^ in the mid 90's. Intensive theoretical 
studies of pairing interactions in nanometric particles followed the experimental work.-"— 
These studies focused mainly on the influence of size quantization of the electronic levels 
on the paired state. The systems considered were relatively large (radius larger than 3nm) 
containing more than 10'^ — 10^ valence electrons, and having irregular shapes, such that 
the only symmetry of the problem is time reversal symmetry. Smaller irregular grains are 
expected to be unpaired, since their average energy-level spacing exceeds the bulk pairing 
gap, breaking the Anderson criterion.-ii^ 

More recently, Kresin and Ovchinnikov^"— investigated the possibility of pairing in 
smaller metallic atomic clusters containing tens of atoms. Such clusters are known to ex- 
hibit a shell structure in their electronic spectrum.—"— This energy structure arises from 
the approximate spherical symmetry of the effective potential felt by the valence electrons 
in the nanoparticle. The near degeneracy of these shells may be considered as a narrow 
peak in the electronic density of states (DOS) which enhances both the energy gap and the 
critical temperature compared to the bulk material. Kresin and Ovchinnikov predicted a 
large enhancement in both parameters (Tc between lOOK to 200K and an energy gap of tens 
of meV) for specific examples of individual aluminum and gallium clusters containing a few 
tens of atoms. Similar results were qualitatively predicted by Friedel in 1992.— The effect of 
shell structure in larger spherical nanoparticles, together with modifications in the effective 
interaction due to alteration of the electrons wave functions, as well as a nonuniform gap 
parameter, were considered in a recent study.— A large and strongly size dependent energy 
gap and critical temperature were predicted for these particles. The effect of shell structure 
on pairing was also discussed in the context of fermionic atoms in harmonic atomic traps.— 

A high DOS at the highest occupied shell (HOS) of the atomic cluster (or nanoparticle) 
is maintained as long as the particle retains spherical symmetry. However, clusters (and 
nanoparticles) with a non-full HOS undergo a Jahn- Teller deformation^S^ which lifts the 
degeneracy of the HOS and reduces the high DOS. Therefore, Kresin and Ovchinnikov^"— 
considered clusters in which the deformation is minimal {i.e. either clusters with completely 
full electronic shells, or with an almost full outer shell) and in which the energy difference 
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between the HOS and the lowest unoccupied shell (LUS) is small, as the most favorable 
scenario for the observation of pairing correlations. 

The energy gap and the critical temperature measured in "large" irregularly shaped 
nanoparticles of various superconductor metals are larger than the corresponding values of 
the bulk material. Evidence for this increase is found in measurements preformed with 
single aluminum particles.-*^ Additional support for this behavior comes from past work 
done with several types of superconducting thin films.—""— These films were composed of 
superconducting grains embedded on top of an oxidized substrate, and separated by the 
dielectric barriers formed by the oxidized substrate. Evidence for a pairing energy gap was 
found in tin grains as small as 2.5nm.— 

The increase in both energy gap and critical temperature was usually attributed to a larger 
effective phonon-mediated electron-electron interaction, due to "soft" surface phonons of the 
individual particles.— i^"— There are also some indications for enhancement of the electron- 
phonon coupling in some molecular devices.—"— However, Kresin and Ovchinnikovi^"— as 
well as Friedel^ assumed a constant phonon-mediated pairing interaction, which was taken 
to be equal to the corresponding values in the bulk material. 

The increase in the energy gap in larger irregularly shaped nanoparticles was also ex- 
plained as being due to finite size modifications in the bulk electronic DOS^i^, and to 
modifications in the bulk electronic wave functions which induce alteration in the effective 
interaction between the electrons.— However, possible changes in the underlying electron- 
phonon interaction, which may be important in smaller nanoparticles, were ignored in these 
works. Furthermore, a semi-classical approximation, together with an assumption of a small 
ratio between the mean level spacing and the bulk gap, were used in these works. In the 
present work, we consider much smaller nanoparticles (containing only hundreds of atoms) 
for which the semi-classical approximation may not be appropriate, and the mean level 
spacing even within the HOS is larger than the bulk energy gap. Therefore, in what follows 
we propose a somewhat different approach which, although limited to simple geometries, is 
more appropriate for the description of pairing in nanoparticles with few hundreds of atoms. 

In this work we present a calculation of the energy gap in nanoparticles at zero tem- 
perature, taking into account the finite size effects on both electrons and phonons. The 
nanoparticles are assumed to be isolated both electrically and mechanically from the outer 
environment. We treat spherical or nearly spherical particles. For this type of particles 
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we explicitly evaluate the electron-phonon interaction, the resulting pairing interaction, and 
the Coulomb interaction. Thus, we are able to predict deviations from the results obtained 
using the average bulk values of these interactions. 

Non-zero temperature and the transition between the paired and the unpaired state 
as temperature is varied are not addressed in this paper. We note that the transition is 
expected to be a "smeared" unlike the one found in an infinite system. Furthermore, the 
usual BCS ratio between the energy gap and critical temperature should not be used in 
order to calculate the characteristic temperature of the transition. 

Three types of materials are considered as examples - aluminum, zinc, and potassium. 
Aluminum is a typical low-temperature, weak-coupling superconductor with the second 
highest critical temperature (in the bulk) of the five superconducting metals - Al, Ga, 
In, Zn, and Cd, that exhibit electronic shell structure in their atomic clusters and ultra- 
small nanoparticles.— Furthermore, aluminum nanoparticles (albeit larger than the ones 
we are considering) are used in experiments examining the properties of superconducting 
nanoparticles. Although potassium is a non-superconducting metal in the bulk its atomic 
clusters and nanoparticles exhibit a pronounced electronic shell structure.— Therefore, the 
effect of the shell structure, together with modifications in the phonon-mediated electron- 
electron interaction, may enable the formation of a paired electronic ground-state. Finally, 
we expect zinc particles to be less susceptible to Janh- Teller deformations than the other 
four superconducting metals mentioned above, due to the larger shear modulus of zinc. 

Pairing was proposed as a mechanism that may explain the apparent odd-even staggering 
in the ionization potential spectra (meaning that the ionization potential of a cluster with an 
even number of atoms is larger than that of the neighboring odd clusters) and in the abun- 
dance spectra of alkali clusters containing up to a few tens of atoms. Pairing interaction 
was also proposed as an explanation for some of the discrepancies between the calculated 
ionization spectrum (obtained by using independent particle models), and the measured 
spectra of sodium clusters containing tens of atoms.— However, it was found that, at least 
when considering these extremely small metal clusters (with up to about 100 atoms), the 
strong fluctuations in the ionization potential could be explained by the deformation of open 
shell clusters^i^ and by the effect of spin degeneracy in the deformed clusters^^ (where the 
only remaining degeneracy is the double degeneracy due to the spin) without resorting to the 
effects of pairing interaction. Therefore, we avoid the regime of very small atomic clusters 
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and concentrate on nanoparticles containing hundreds of atoms. 
The general scheme of our calculation is as follows: 

1. Using a simple and crude model, the single-electron and single-phonon states and 
spectra are evaluated. 

2. The electro n-phonon interaction is calculated assuming a static Thomas- Fermi 
screened interaction. 

3. The phonon- mediated pairing interaction is evaluated using either the Frohlich trans- 
formation or the similarity renormalization. 

4. The splitting of the electronic degenerate levels caused by deviations from spherical 
symmetry of the particles is estimated. 

5. As a reference, the modification of the electronic spectrum due to the pairing interac- 
tion is calculated in spherical particles. 

6. The reduction in the effects of pairing on the low energy excitations of the electronic 
spectrum due to small deformations is evaluated. 

7. Finally, the reduction of the pairing interaction due to Coulomb repulsion is also taken 
into account. 

We note that the effective interaction is initially calculated by applying the Frohlich 
transformation^ to the system's Hamiltonian. This method is easy to apply and it describes 
correctly some of the features of the pairing interaction. However, it produces unrealistic 
results when applied to a system with non-degenerate electrons. It therefore can be used 
only for the spherical particles. In order to account for the effect of deformations, we use the 
similarity renormalization method.-^ This method is known to give accurate results for the 
transition temperature and the energy gap in strong coupling bulk material superconductors 
within the framework of a Hamiltonian theory.— 

The pairing Hamiltonian of the deformed particles involves varying coupling coefficients 
and non-uniform energy level distribution. We are unaware of an exact solution to such 
a Hamiltonian. Therefore, we use the simple BCS grand-canonical approximation to in- 
vestigate this Hamiltonian, although, in principal, a canonical ensemble approach, such as 
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fixed-N projection of the BCS state^*^'^'^, is more appropriate. Several authora^i^"— have 
pointed out that finite level spacing should modify the BCS solution, leading to a smaller 
lowest pairing excitation energy compared to the one obtained by the BCS approximation. 
This reduction is found to be approximately equal to the level spacing of the single-electron 
spectrum as long as the relevant level spacing (the level spacing within the HOS in our 
model) is smaller than the lowest pairing excitation energy. In our model we find that the 
lowest excitation energy is several times larger than the HOS level spacing, and therefore 
this finite-level-spacing correction can be neglected as long as one is merely interested in a 
rough estimate of the modifications to the lowest excitation energy. Therefore, for our pur- 
poses the simplicity of the BCS approximation makes it an adequate choice for calculating 
the excitation energy (see some additional comments on this subject in section IVIBp . Nev- 
ertheless, a more accurate description of the pairing correlations should involve a canonical 
ensemble treatment. 

The structure of the paper is as follows. In section [TTl we describe the single particle states 
of both electrons and phonons in spherical particles, as well as the interaction between them. 
The derivation of the effective electron-electron interaction and the resulting approximate 
model pairing Hamiltonian is given in section IIIII The deformation of the spherical particles 
is estimated in section HVl while the effect of Coulomb repulsion is discussed in section |Vl 
The modification of the free electrons spectrum is addressed in section IVll Numerical results 
concerning both the degenerate and the non-degenerate scenarios are presented in section 
IVIIi Some measurable quantities, which may be used to identify pairing in such small 
nanoparticles, are briefiy discussed in section IVIIIi Our main conclusions are summarized 
in section HXl 

II. ELECTRONS, PHONONS, AND THE INTERACTION BETWEEN THEM 
A. Single electron states 

The delocalized electrons of a nanoparticle are modeled as free particles within an infinite 
spherical potential well. Alternatively, one may use a spherically symmetric harmonic po- 
tential. The former effective potential corresponds to complete screening of the electrostatic 
potential of the particle's positive ionic background, while the latter reflects a complete lack 
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of screening. As we comment later, it is reasonable to assume that the screening picture of 
the bulk material essentially persists even in nanoparticles that contain hundreds of atoms. 
Thus, we may expect that the electrostatic interaction is almost completely screened. Fur- 
thermore, the major features of the shell structure observed in the abundance and ionization 
spectra of atomic clusters, reflect the filling of major energy shells that are characterized 
by an angular quantum number / and a radial quantum number n.—^ This fits well with 
the spherical box approximation, although a deformed 3d oscillator gives an even better 
fit.— 1^ The realistic effective potential lies somewhere between these two extremes, and 
includes deviations from spherical symmetry.— Note, however, that the 3d harmonic po- 
tential introduces an artificially high degeneracy of the energy levels, which would lead to 
an enhancement of the pairing effects we are exploring. This additional degeneracy is lifted 
even by a small amount of screening. Thus, the infinite well approximation captures better 
the physics of the true potential. 

The energy levels of electrons in a spherical box are given by 



and their corresponding wave functions are 



where / and m are the usual angular momentum quantum numbers, with / > and — / < 
m < I. Also, ji are the spherical Bessel functions, is the nth zero of the function ji 
divided by the radius R of the sphere, m* is the effective mass of the electrons, and Yim are 
the spherical harmonics. We assume that the effective mass of the electrons remains the 
same as in the bulk metal.— The normalization constant Bin is given by 

Bin = ^{3i+i{kinR))-\ (3) 

We note that, including spin degeneracy, each energy level may contain up to 4/ + 2 electrons. 
The electron field operator is given by 






(4) 



Imna 



and the second quantized Hamiltonian takes the usual form 




(5) 



Imna 
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where a is the spin index. 
B. Single phonon states 

Acoustic phonons are described as elastic disturbances in the positive ionic background 
of the nanoparticles, which we model by using linear elasticity theory. Since the typical 
nanoparticles we are dealing with contain few hundreds of atoms, the use of an elastic 
description could be questioned. However, the applicability of elasticity theory for nanopar- 
ticles is supported by a reasonable description of the low frequency acoustic modes for 
nanometer size particles.—"— Also, phonon related phenomena in nanomechanical systems 
are often modeled by using an elastic model for the description of the vibrational degrees of 
freedom.—"™ Therefore, although in principle one may solve the discrete equations of mo- 
tion for the few hundred atoms, we prefer the simple elastic model that enables an analytic 
expression for the effective interaction between the electrons. 

The small size of particles under consideration leads inevitably to effective elastic con- 
stants that deviate from the bulk material values. Generally, nanometer size systems tend 
to exhibit smaller elastic coefficients compared to the bulk. However, ab-inito density func- 
tional theory calculations of the effective modulus of dilation of small silicon, tin, and lead 
atomic clusters and the effective Young modulus of ultra-thin Si nanowires^*^, indicate 
that the softening of the elastic constants is quite restricted. For example, the calculated 
reduction of the effective modulus is about a factor of two for silicon clusters containing 
as few as 15 atoms, and about 25% for tin and lead particles with the same number of 
atoms. The calculated variation in the Young modulus of silicon wires is even smaller. We 
also note that experimental studies have found that the Young modulus of gold nanorods 
(with diameters varying between 10nm-20nm) is either lowered by 20%-30%^^'^^ or remains 
essentially unchanged^ compared to the bulk material. We therefore use both the Lame 
constants of the bulk material, as well as 25% smaller Lame constants for comparison. 

We solve the linear elasticity equation of motion under stress-free boundary conditions 
imposed on the surface of the sphere. The relevant results of the elastic solution are given 
in appendix [X] based on the treatment by Eringen.™ The normal modes of vibration can 
be divided into two types: spheroidal modes containing both longitudinal and transverse 
components, which are coupled through the boundary conditions, and torsional modes con- 
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sisting solely of transverse components. Similar to the electronic wave functions and energy 
levels, the vibrational modes are characterized by two angular momentum quantum numbers 
/ and m, and by a radial quantum number n. 

As will be explained in the following section, we are interested only in the longitudinal 
component of the spheroidal modes. This component can be written as 

P 

where cj) is the solution of the scalar Helmholtz equation 

(V2+p')0 = O, (7) 
which for stress-free boundary conditions takes the form 

(t)ljnn = Anjl {PlnT) Yim{9 , 0), (8) 

where pin = w^n/cio (cio is the bulk longitudinal sound velocity of the material). The bound- 
ary conditions impose a discrete spectrum, where the eigenfrequencies are characterized by 
the angular quantum number / and by the radial quantum number n. Furthermore, the ra- 
tio between the amplitude Ain of the longitudinal component and the amplitude Cin of the 
transverse component are uniquely determined by the boundary conditions'^ (see appendix 
We note that a continuum approach results in an infinite number of normal modes. In 
order to account for the atomistic discreteness of the nanoparticles, we take into account 
only the lowest 3Na modes of the entire spectrum, where A^^ is the number of atoms in 
the nanoparticle. We define the Debye energy of a nanoparticle as the energy of the most 
energetic phonon mode. 

In order to quantize the vibrational modes we write the spheroidal displacement field as 




where p is the material density of the cluster and 6^ (6) are the phonon creation (annihila- 
tion) operators. The second quantized Hamiltonian is obtained by inserting the quantized 
displacement field into the linear elastic Hamiltonian, thus obtaining 

Hph = y^Lnblmn + ^ j • (lO) 

Imn ^ 



9 



The amplitude of the longitudinal component 
tion condition imposed by the quantization procedure 



and by the ratio Ain/Cin- 



IV 



sph 
Imn 



Ain is uniquely determined by the normaliza- 
ocedure 

'd\ = l, (11) 



C. Electron-phonon interaction 

We assume a screened interaction between the electrons and the density disturbance in 
the ionic background, i.e. the phonons. This interaction is added to the free Hamiltonian 
composed of and Hph. We take into account only the static electronic screening when 
considering the interaction between the electrons and the phonons.—"— This approximation 
is reasonable since the maximal frequency of the phonons is two orders of magnitude smaller 
than the Fermi energy of the electrons. In order to keep the calculation as simple as possible 
we assume static electronic screening also when considering the direct Coulomb interaction 
between the electrons (see section |V]). We apply a unitary transformation to the basic 
Hamiltonian {He + Hp + He-p) in order to account for the contribution of the phonons (or 
the phonon screening) to the effective interaction between the electrons, as in the original 
work of Frohlich.— 

We utilize the Thomas- Fermi approximation in order to describe the static screening. This 
approximation should still be reasonably valid for particles containing hundreds of atoms, 
since the typical bulk Thomas- Fermi screening length is an order of magnitude smaller than 
the diameter of such particles. Therefore, the screened electrostatic interaction Hamiltonian 
between the electrons and the phonons takes the form 

He-p= peiri)vTF{\ri-r2\)5pi{r2)d^rid^r2, (12) 

Jv Jv 

where 

VTF{r) = 

r 

oo I 

= AnkTF J2 i^TFr<) h {kTFr>) J] 1/™ (^i0i) Y^l (^2^2) , (13) 

1=0 m=-l 

Pe is the electronic local density, 6pi is the local variation in the ionic background charge 
density due to the presence of the phonons, /ctf is the Thomas- Fermi wave number, ii{kr) 
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and ki{kr) are the modified spherical Bessel functions, r> = ri and r< = r2 when ri > r2, 
and vice versa. 

Considering small deformations, the relative change in an infinitesimal volume element 
is given (to first order in the displacement field) by V-uJ^ Therefore, to first order in u, 
the effect of the transverse component of the spheroidal phonons on 6pi vanishes. Thus, 
we can take into account only the longitudinal component of the spheroidal modes when 
considering the interaction with the electrons. 

The interaction Hamiltonian is therefore given by 

i/e-p= / / e2Zno^t(ri)^(ri)i;TF(|ri-r2|)V-u^P^(r2)rfVirfV2, (14) 
Jv Jv 

where Z is the number of valence electrons per atom, and uq is the atomic density of the 
cluster, which is taken to be equal to its bulk value. Using Eqs. (I2])^(I3]), (in])^(in]), and f lTT]) 
we obtain 

He-p = MLMNcl_^rninia'^l2m2n2abl3m3n3 + C.C, (15) 

LMNa- 

where L = {li, 12,13}, M = {1711,1712,7113}, and N = {^i, 7^2,^3}, and where the angular 
momentum quantum numbers satisfy 

k > 0, i = 1,2,3 (16) 
— /j < rrii < li, i = 1,2,3 (17) 
1/1-/2! < h<h + l2 (18) 
h + h + h = even integer (19) 
^3 = i7ii — ni2. (20) 
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The coupling coefficient Mlmn is given by 

Mlmn = HLNc{li,l2,h\rnx,-m2,mi- m2){-l)'^^ 

= CLNRLNQLc{ll,l2,h;mi, -1712,7711 - 1712) (-1)'^^ (21) 



R 



Rln = I drirji^ (ki^mn) ji^ {ki^nin) 







ri 



ki3 (kTFri) / rfrsrsJig (p«3„3r2) % {,kTFr2) 







+ii3 {kTpri) I dr2r2ji3 {,Phn-,r2) {kTFr2) 



(23) 



where c (/i, /2, ^3; 0, 0, 0) and c{li,l2,h',77ii, —m2,mi — 1712) are Clebch-Gordan coefficients 
(adopting the notation of Rose-^). The detailed derivation of M^mn is given in appendix 
IBI We can summarize the total Hamiltonian as 

H = He + Hph + He-p = Hq + He-p. (25) 



III. EFFECTIVE ELECTRON-ELECTRON INTERACTION 
A. The Frohlich transformation 

We apply the unitary Frohlich transformation'^^ to the Hamiltonian (l25!l . 

= e^'He' = H, + H,^p+[H,, s] + ^ [{H,.p + [Ho, s\) , s] + ^ [H,.p, s]+0 [MI^j^] , (26) 
where s is an anti-hermitian operator defined by 

LMNa 

The coefficients aLMN are given by 

C^LMN — ^ — T ) \^°) 

and are chosen so as to eliminate the electron-phonon interaction up to second order in 
Mlmn- Specifically, the particular choice of a ensures that ifg-p + [Hq, s] = 0, and thus the 
lowest non-vanishing term in the transformed Hamiltonian is | [-f^e-p, s] . 



12 



The interaction term contains a sum of operator products, most of which still involve the 
phonon operators. However, it also contains a purely electronic term, which represents the 
effective electron-electron interaction resulting from the original electron-phonon interaction 
(up to second order in Mlmn)- Calculating explicitly this interaction term we obtain the 
following effective electron-electron interaction 



^^" = E E E 



l-i ioi'^'Tii momi —mo / ,' , , ' ' ''o"'o 



/ / 



where the primed indexes stem from the operator s. The summation over is restricted 
only by the finite number of modes in each branch of the phonon spectrum (defined by the a 
specific value of /s), due to the finite number of normal modes supported by the nanoparticle. 
We note that the numerator in fl29l) is not always positive, unlike the numerator appearing 
in the similar result obtained for translationally invariant system.— 

Following Cooper's argument?^, the electrons in the paired state are paired in a manner 
that "uses" the attractive part of H^-e in the "most efficient way". First, we note that 
the difference between adjacent electronic energy levels is typically larger than the maximal 
typical phonon energy by an order of magnitude or more. Therefore, we can expect to 
obtain a maximal attractive interaction between electrons belonging to the same energy 
level. Thus, we neglect inter-level phonon-mediated interaction. Next, we assume, as in 
the ordinary BCS theory, that the electrons are paired in a singlet state implying a' = —a. 
Finally, due to the properties of the Clebch-Gordan coefficients appearing in (Elj), we find 
that, in order to ensure the negativity of the intra- level interaction, we must assume that the 
Cooper pairs are composed of electrons with opposite z component of angular momentum. 
This restriction means that 

mi — 1712 + fT^2 — (30) 
1712 = (31) 
mi = —m . (32) 

Eqs. (l30!) - (l32l) ensure that the numerator appearing in (129|) is always positive for time- 
reversed pairs of electrons of the form {m 4-}- Since inter-level interaction is ignored, 
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all fully occupied levels, or levels with one electron or one hole, are considered to be inert 
and the phonon-mediated interaction acts only within the HOS. This approximation holds 
as long as the temperature of the nanoparticle is much smaller than the level spacing around 
the Fermi level. Since this level spacing is typically of the order of O.leV for the particles 
we are considering, the approximation is reasonable up to a temperature of a few hundred 
Kelvins. Therefore, we are left with the effective Hamiltonian 

|2 



mm' hn3 



I Trnnn-j | 



1 2 



He-e = XI XI ^ ^' -m ,m - m) (-1)"^ ^L'naim'n-a^'-"^"-'^^'"'"'^' 



(33) 

where / and n characterize the highest occupied electronic level. 

The structure of the effective Hamiltonian in f p3|) may have been anticipated in advance, 
since the coupling between time-reversed pairs is a common characteristic of the pairing 
phenomenon.— 1^ Our numerical calculations for the three materials under consideration, 
show that the coefficients aiMN in (l28ll are always much smaller than 1, when one considers 
the interaction within a degenerate HOS. Therefore, neglecting higher-order terms in ( 126|1 
is justified. 



B. The similarity renormalization 

In this section we introduce an alternative derivation of the interaction Hamiltonian. The 
above calculations work well for exactly degenerated HOS. However, in the next section we 
want to study the effects of the splitting of the initially degenerate electron energy levels 
due to static deformations of the spherical particles. We expect the splitting to reduce 
the effective electron-electron interaction compared to the spherical system. The effects of 
splitting cannot be accounted for correctly in the framework of the Frohlich transformation. 
Using f l29|) and introducing the energy splitting due to deformations, one obtains an artificial 
enhancement (and even divergence) of the effective interaction, when the energy splitting 
e^ini — ^i2n2 is of order hwi^ns- Note that even if no actual divergences are encountered, the 
second order expansion ( l26ll is invalid if the denominators in ( l28ll are smaller than the values 
of Mlmn- 

Another deficiency of the Frohlich interaction fl29|) lies in the fact that it includes terms 
that represent coupling of electronic states by a virtual phonon whose energy is smaller than 
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the energy separating the electronic states. In the usual BCS treatment these (repulsive) 
terms are avoided, since the interaction is assumed to be constant in fc-space, with an 
artificial cutoff at the Debye energy of the material.— Lastly, using a Frohlich type interaction 
within the BCS formalism results in a large overestimate of the size of the energy gap and 
critical temperature of bulk superconductors.— 

An alternative derivation of the effective interaction relies on the application of the sim- 
ilarity renormalization method to the initial Hamiltonian fl25|) . Although its application is 
more complex, it avoids the appearance of vanishing energy denominators. Furthermore, 
the resulting interaction is reduced between electrons with different energies, and it is auto- 
matically cut-off at the right energy scale. Also, the obtained effective interaction is always 
attractive. As shown by Mielke^i^ this approach yields the correct critical temperature 
and energy gap for strong-coupling bulk superconductors within the framework of the BCS 
model. 

The derivation of the effective interaction is given in appendix [Cl where we follow the 
treatment by Mielke.~ We obtain the following effective electron-electron interaction be- 
tween electrons belonging to the same energy shell 

I 1 2 

2 ]\/[ , 

Gram' = ^ TZ _ ^^T^ © (^^3^3 " l^mhna " ^m'hns I) ' (^4) 

where G is the Heaviside function and Smhna are the single-electron energies in the HOS 
of the spherical or deformed nanoparticles. The Heaviside function in fl34|) ensures that a 
phonon can mediate interaction between two electron states only if its energy is larger than 
the energy separation between the two states. We note that ( l34l) and ( !33l) coincide if the 
deformations are neglected. 

Corrections to the wave functions of both electrons and phonons due to deviation from 
spherical symmetry are neglected in flM|) . This is justified because it affects the numerator 
in (!34l) only to second order in the perturbation. The effect of deformation on the phonon 
spectrum is minor, as deformations are small compared to particles' radii (see section IVlip . 

Our result differs from Mielke's interaction in the cutoff function appearing in G^^i 
( l34|) . The source of the difference between our cutoff function and Mielke's cutoff function 
is explained in appendix O However, our interaction and Mielke's coincide for the physical 
scenario considered by Mielke. This scenario consists of an effective interaction mediated 
by non- dispersive Einstein phonons, whose frequency is much larger than the electron en- 
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ergies. Our result is similar to the one obtained by Hiibsch and Becker— using a different 
perturbative renormalization scheme. 

Finally, we note that the electron energies are renormalized by the electron-phonon in- 
teraction, and that the effective interaction should depend on the renormalized electron 
energies.— i^i^ The renormalization is especially important in describing properties of strong- 
coupling bulk superconductors.— 1^1^ This is another deficiency of the Frohlich interaction, 
which depends on the non-renormalized energies. By contrast, the energies appearing in Eq. 

are in fact the renormalized energies. However, renormalization of the single-electron 
energies is not important in our model since the electron-phonon interaction cannot lift the 
degeneracy of the electron shells in the spherical particles. The contribution of renormaliza- 
tion to the splitting of the HOS in the deformed particles is small, since it is of the order 
of a]^2 [M'^/hu), where ^ parameterizes the deviation from spherical symmetry 

[the definition of a]^2 is given in Eq. (H2|) in the following section]. Therefore, we neglect 
the effect of renormalization on the electron energies and instead use the non-renormalized 
electron energies. In appendix O we show that the renormalization of the phonon energies 
due to the electron-phonon interaction can also be neglected. 

The energy shift due to renormalization does depend on the shell quantum numbers I and 
n and the electron filling. Therefore, in a more detailed treatment that takes into account 
inter-shell effects (like the one carried out by Kresin and Ovchinnikov^"— ) renormalization 
effects should be considered, especially when dealing with nearly degenerate HOS and LUS. 
In appendix O we give the expression ( 10231) for renormalization of electron energies due to 
electron-phonon interaction within the HOS, which is responsible for most of the HOS shift. 
The expressions for the renormalization of electron energies and phonon energies are similar 
to the ones obtained for the bulk system by Mielke.— 



IV. DEVIATION FROM SPHERICAL SYMMETRY 

Deviation from spherical shape leads to a splitting of the degenerate electron levels and 
to a decrease in the total energy in the free electrons of the nanoparticle. On the other hand, 
deviation from a spherically symmetric shape leads to an increase in the elastic energy of 
the nanoparticles. The magnitude of the deviation is determined by the balance between 
the increase in the elastic energy and the decrease in the electronic energy. A similar type of 
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calculation was used by Kresin and Ovchinnikov^ in order to estimate shape oscillations of 
aluminum atomic cluster (with 14 atoms), in which an electron from the full HOS is raised 
to the LUS. 

The general expression for the elastic energy density is^ 

2 



1 , ,9 [1 / dui dun\ 



(35) 

J L \ J * / J 

where X,fi are the Lame coefficients, and u is the displacement field. For simplicity, we 
assume that the local density in the nanoparticles remains the same as in the bulk material 
regardless of the shape of the particles. Therefore, we need to consider a distortion of the 
spherical shape for which 

V ■ u = 0. (36) 

Choosing such a displacement field eliminates the first term in ( 135|) . We also assume an 
axially symmetric distortion of the surface of the spherical particles. This distortion can be 
parameterized as follows 



R 



cos( 



(37) 



where Pi{cos6) are the Legendre polynomials. A displacement field that produces such a 
distortion while not changing the local density is given by^ 

u = ^ V {r'Pi (cos 6)) . (38) 

I 

We are interested only in a rough estimate of the magnitude of the deformation. There- 
fore, we do not consider the general displacement field f l55]) . Instead we only examine 
separately each value of I in (1381) and determine which multipolarity enables the largest 
decrease in the total energy compared to the spherical particles. 

For a certain multipolarity I, the first order correction to the energy of a given electronic 
state (characterized by the quantum numbers L, M, and N) is^ 

S^LMN = -'2eLNai < LMN\Pi{cose)\LMN > 

= -2eLNaic{l,L,L]0,0,0)c{l,L,L]0,M,M) . (39) 

Only even values of < / < 2L yield a nonzero 5e^^j^,ij^. For a full shell the sum ^j,^ Sclmn 
vanishes, and therefore we need to consider only non-full HOS. 
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For a given / the additional elastic energy is 

5E[i, = / £reii, = ciTia^R^fx, (40) 



where q is a /-dependent constant that hes between 2 for / = 2 and 4 for large /. We assume 
the deformation to be small and thus integrate the radial functions between zero and the 
original radius of the sphere. Consistently, we find the resulting deformation to be small, 
with a maximal aspect-ratio of 1.13 and an average aspect-ratio of about 1.05. 

AEi=J2S^iMN + SKla, (41) 
M 

where the summation in fl4ip is over all the occupied states in the HOS. The states are filled 
from |M I = and upwards until the highest relevant value of |M|. 

We obtain the largest decrease in the energy of the nanoparticles Ai?™"^ and the cor- 
responding amphtude by differentiating AEi with respect to ai and equating the 
derivative to zero. The resulting expression for the amplitude is 

,max ^LN 



a, 



c (/, L, L; 0, 0, 0) ^ c (/, L, L; 0, M, M) . (42) 



The maximal value of is obtained for half filling of the HOS. 

The largest |Ai?™°^| is obtained for / = 2 (quadrupole deformation). |A£';™2^| is larger 
by a factor of 2 to 10 than the second largest |Aii^["°^| (for the relevant values of L), and 
larger than the sum of all other lAE'^^^gil with their corresponding optimal a^"^. 

We ignore the effect of surface tension which tends to decrease the deviation from the 
spherical shape. This effect is probably important because a large fraction of the atoms 
in the nanoparticles we study reside on the surface of the particles. Taking into account 
only the bulk elastic energy and ignoring the increase in surface tension we overestimate the 
size of the deformation and of the energy splitting. We also note that the deviation from 
spherical symmetry obtained here is not much larger than the effect of surface roughness 
due to the discrete atoms. However, surface roughness corresponds to a large / deformation 
which results in a small energy correction Se^j^j^j^ compared to the quadrupole correction. 



V. EFFECT OF COULOMB REPULSION 



We use the same type of screened Thomas-Fermi potential to describe the repulsive 
Coulomb interaction between the electrons. The Coulomb part of the Hamiltonian can be 
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written generally as 



^co — ^ ^L°MN'44m4n4cr42m2n2a' '^hminirr' '^hrnsn^a, (43) 



LMNaa' 



where L = {li,l2,h,U}, M = {mi, m2, ms, 7714}, = {ni, ?t,2, 713, 774}. The Coulomb inter- 
action matrix element is denoted by Ml°^^j^ 

(44) 

We focus on the Coulomb matrix elements that are relevant for the interaction between 
time-reversed electron pairs within the same energy shell. We insert the expressions of the 
electronic wave functions (|2]) and of the screened Thomas-Fermi interaction (fT3|) into the 
Coulomb matrix element (HI]), and obtain 

R^JT+i {huR) Jo 

rri rR 

k^' {kxpri) I dr2rljf {kinr2) i^' {kTFr2) + «l' (kxpri) / dr2rljf {kinr2) k^i {kTFr2 

Jo Jri 

j dn.Y^i, [n,) Yi„, {n,) y^,,,, {n,) j dQ2Y;_^, m Yi_„, {^2) r;,^,, (^2) ■ 

The angular integrals yield 

d^iY^^iYimYj^i j^^i = 

_^±i__c(/, /, L'; 0, 0, 0)c(/, /, L'; m, -m', -M') (46) 

(-I)-'"' ^Jii^=c(/, /, L'; 0, 0, 0)c(/, /, L'; m, -m', -M'), (47) 

where c{l, /, L ; 0, 0, 0) = unless L is an even integer, which means that (-1)^ = 1. Using 
Eqs. din]) and f H7)) . one obtains the Coulomb part of the Hamiltonian 

H,o=Ni^Yl E J2^^nL'Q^L'[c{l,l,L';m,-m',m-m)] (-1)"' 

m,m'=—lL'=0 



(45) 
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where, 



= Jfn'. (49) 



rR 

RinL' = / drirljf (kinri 
Jo 



R 



+iL' {f^Tpri) / dr2r^ji {hnT^) k^> {kTFr2 



(50) 



{21 + 1) 
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We cannot eliminate the Coulomb interaction between electrons within a degenerate 
or nearly-degenerate HOS by adding an additional anti-hermitian generator to the Frohlich 
generator (1271) . since such a generator contains vanishing or very small energy denominators. 
However, if we transform the Hamiltonian (l25l) (where we replace H^-p with He-p + Hco) 
using the Frohlich generator fl27|) . we obtain the expansion 



e^'He' = Ho + H,,+ ^ [H,_p, S] + [H,o, S] + .... (52) 

The leading correction to the electron-electron interaction due to Coulomb interaction is 
just Hco- The commutator [Hco, S] does not renormalize the electron-electron interaction. 
Therefore, the next correction is of higher order in the commutator expansion, [[Hco, S] , S], 
i.e. of order M'^M^°. We neglect these higher order terms and take into account only the 
lowest order correction, Hco. A similar result is obtained if the similarity renormalization 
is applied to an initial Hamiltonian containing the Coulomb term and the generator of the 
transformation is constructed to eliminate only the electron-phonon interaction (for more 
details see appendix [C]) . Thus, in all cases considered we simply add the screened Coulomb 
term to the phonon-mediated interaction. 

The effect of Coulomb interaction on pairing correlations in the bulk is significantly 
reduced due to renormalization of the average Coulomb interaction constant.— The renor- 
malization in the bulk stems from the large difference in the energy scale in which the 
phonon-mediated interaction acts (up to the Debye energy) and the one in which the screened 
Coulomb interaction acts (up to the plasma energy of the material which is approximately 
the Fermi energy) , and as well as from the fact that one can reasonably well approximate the 
Coulomb interaction, over the entire energy range, by a single average constant. However, 
in our system, only relatively narrow energy shells contribute to the renormalization of the 
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Coulomb interaction. Furthermore, the inter-shell Coulomb interaction matrix elements are 
much smaller than the intra-shell matrix elements. Therefore, they should not contribute 
much to the renormalization of the intra-shell Coulomb interaction. We therefore claim 
that the relevant energy scale for the action of the Coulomb interaction in our system is 
the width of the split HOS and not the Fermi energy of the particles. Accordingly, we can 
take into account only the Coulomb interaction within the HOS, and the renormalization of 
the average Coulomb interaction is expected to be smaller than the renormalization in the 
bulk. A similar situation occurs in superconducting Ceo, where the Coulomb pseudopoten- 
tial constant is not much smaller than the non-renormalized average Coulomb interaction 
constant. '^^ By using the detailed Coulomb interaction matrix elements and not an average 
interaction, we avoid the necessity of estimating the renormalization within the HOS. 



VI. MODIFICATION OF THE ELECTRONIC SPECTRUM DUE TO PAIRING 

We use two types of approximations to study the manner in which the highest occu- 
pied electronic level is split due to the pairing interaction at zero temperature. In the 
first approximation, which is applicable only to a degenerate HOS, we average the effective 
electron-electron interaction over the entire HOS, and obtain a single electron-electron cou- 
pling constant. This model (often referred to as the "seniority model") has the advantage 
that it is analytically solvable for a fixed number of electrons in the HOS.— The second 
approximation consists of using the BCS grand-canonical approximation, in which we fix 
the average number of electrons to be equal to the true number of electrons in the HOS. 
In this approximation we can relax the requirement of a single coupling coefficient and of a 
degenerate HOS. 

We note that the seniority model was recently used^^ to evaluate several properties (such 
as temperature-dependent specific heat and magnetic susceptibility) of spherical nanoparti- 
cles, in which a constant pairing interaction extrapolated from the bulk material was assumed 
to act within a completely degenerate HOS. The effects of a uniform magnetic field were 
also considered. 
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A. The seniority model 



We replace the coupling coefficients of the Hamiltonian (l33i) with a single average coupling 
constant 



G 



1 I ^"""3 1 2 2 

TTiyiyi^^ c{l,l,l3;m,-m' ,m-m') . (53) 

' mm' hns 



(2^ + ir±:ie ^"3^3 



The pairing Hamiltonian ( 133|) is replaced by 



The energy levels of the Hamiltonian flMl) are given by 



E(^,Ar) = ^(A^-5)(2^]- Ar-5 + 2), (55) 



and their degeneracy is given by 



D{S) = { 



1, 5 = 

2Q, S = 1 (56) 



20,1 20! 



S>2, 



{2n-S)\S\ (2r2-5+2)!(5-2)!' 

where is the total number of electrons in the HOS, Q is equal to 21 + 1, and 5* is the 
"seniority number", which counts the number of unpaired electrons in the HOS. In other 
words, S is equal to twice (twice plus one) the number of broken Cooper pairs in the HOS 
if N is even (odd). When the number of electrons in the HOS lies between 2 and 21 + 1, the 
seniority number has the following values 

f 0,2,...,iV AT even 
S = < (57) 
[ 1,3,..., AT iVodd. 

The expression ( l55l) for E{S, N) can be used even if the number of electrons exceeds Q, but 
then should be taken as the number of holes in the HOS.— 

The energy difference between two adjacent levels with seniority numbers 5* and 5* + 2 is 
equal to 

E{S + 2,N)-E{S,N) = -G{Q-S). (58) 

Thus, the energy levels become denser when considering higher values oi S {i. e. higher energy 
levels). We consider the energy difference between the ground-state {S = or S = 1) and 
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the first excited state (S* = 2 or S* = 3) as the 'energy gap'. The energy gap is equal to —Gfl 
if the number of electrons in the HOS is even, and to —G{Q — 1) if the number of electrons is 
odd. Note that, unlike bulk material, the size of the gap varies linearly and not exponentially 
with the magnitude of the coupling constant. This result remains approximately valid even 
when deformations are taken into account. 

B. The BCS model 

The analytical results of the seniority model cannot be used for a non-degenerate HOS. 
Therefore, we employ the BCS approximation in order to analyze the pairing Hamiltonian 
in the non-degenerate HOS. The applicability of the BCS approach was discussed in section 
m Here we add some comments about this subject. 

The BCS approach was found to be reasonably successful in describing the ground- 
state and low-excited states for systems containing tens of interacting fermions, such as 
the nucleus.-*^"— In particular, Braun and von Delft^ conclude (based on literature dealing 
with pairing in the nucleus) that the BCS approximation is adequate in order to describe, 
at least qualitatively, pairing correlations in ultrasmall nanoparticles. Indeed, considering 
the lowest excitation energy of pairing Hamiltonians in the context of nuclear pairing, the 
deviation between the BCS results and the ones obtained by more sophisticated treatments 
or the exact solution was found to be up to two-fold^i^, while the typical differences, both 
in the context of nuclear pairing^i^ and electron-pairing in metallic nanoparticles^"^, are 
usually smaller. We may conclude that following these comments, together with the discus- 
sion in section HI the usage of the BCS approximation is appropriate within the framework 
of our model, as long as we only aim at a rough description of the low excitation energies of 
the pairing spectrum of the nanoparticles. 

We note that, in what follows, we do not assume the usual constant pairing interaction and 
gap parameter as in the usual BCS approximation. Instead we use the detailed interaction 
matrix elements, and solve for the m-dependent gap parameters A^. 

The BCS gap equations for a given HOS, characterized by quantum numbers I and n, 
are given by (see for example Greiner and Maruhn^) 




(59) 
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where Em are the energies of the electrons in the HOS as a function of the angular quantum 
number m, A is the chemical potential, and G^^' is the pairing potential acting between 
the electron pairs {m, —m} and {m', — m'}. The coupling coefficients are given by 

G. 



o TT-rin-i 



2 

c{l,l,l3;m,—m,m — m) , (60) 



if deformations are ignored, and by Eq. if they are taken into account. The set of 2/ + 1 
equations defined by fl59|) is solved together with the equation 

AT = y ( 1 _ , - ^ I . (61) 



Eq. fl6T|) is obtained by setting the expectation value of the number operator at the BCS 
ground-state to be equal to the total number of electrons in the HOS. We solve the BCS gap 
equations only for the case of an even number of electrons in the HOS, since the expectation 
value of the number operator in the BCS ground-state is an even number. In correspondence 
to the definition of the energy gap in the seniority model, we define the energy gap as twice 
the energy of the lowest BCS quasi-particle. The energy of the BCS quasi-particles is given 
by 



= ^JAI + {e^ - Xy. (62) 

Considering only the phonon-mediated interaction in the degenerate case ( l60l) one finds 
that the sum G^^' does not depend on m. Thus, obtained by solving Eq. ( 15^ 
is independent of m (although the G^^' are not equal to each other) and all BCS quasi- 
particles have the same energy. In this particular case, the energy gap obtained by the 
seniority model (for an even number of electrons) is reproduced by the BCS model. The gap 
parameter becomes truly m dependent only for non-average m-dependent interaction 
matrix elements, and non-degenerate HOS. 

By considering deformation together with the Coulomb interaction, we obtain for a large 
portion of the particles an overall repulsive average interaction. In such cases the standard 
mean-field approximation of the gap parameters. Am = A, yields an unpaired ground-state. 
However, we do not apply this approximation and therefore are able, in principal, to find 
non-trivial solutions to the gap equations. In fact, even for an entirely repulsive interaction 
such solutions may exist as long as the interaction is not constant. 
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The scenario of pairing in the presence of a repulsive interaction was first considered by 
Tolmachov^, and since then in numerous variations by many authors.—"— Furthermore, 
Mila and Abrahams^fl have shown that, for a bulk superconductor, a solution to the gap 
equation which is an odd function of k — kp is possible for an arbitrarily strong short-range 
repulsive interaction, and that this solution has a lower energy than the unpaired state. 
For the aluminum and zinc nanoparticles we find that at least a small part of the overall 
interaction matrix elements is in fact negative. 

The energy of the paired ground-state (relative to the new chemical potential A) is 

E^.. = j:(e^-X)(l- V^^./l-^fl^ijgZ. (63) 

V y^;. + (e™ - A) V '» » A J. + (e„.-A) 

and the condensation energy of the nanoparticles is defined as the difference between Epair 
and the energy of the free-electron ground-state. The gap equations ( l59i) may possess sev- 
eral solutions. Of those solutions we choose the one with the lowest condensation energy. 
However, since the different gap parameters may have different signs, it is possible that 
the condensation energy of this solution is positive. Only if the condensation energy of the 
solution with the lowest pairing energy ( l63l) is negative we can claim that the paired state 
is stable and that it represents the true ground-state of the nanoparticle. We find such 
solutions for a large portion of the aluminum and zinc nanoparticles with non-degenerate 
HOS (see section |VllD]) . 

VII. RESULTS 
A. Main results 

We consider aluminum nanoparticles containing 100 to 400 atoms (or 300 to 1200 free 
electrons), zinc nanoparticles containing 200 to 500 atoms (or 400 to 1000 free electrons) and 
potassium nanoparticles with 100 to 500 atoms. Of the three metals considered, our results 
indicate that aluminum is the best candidate for observing pairing effects in ultra-small 
nanoparticles containing a few hundreds of atoms. 

The average electron-electron interaction in aluminum nanoparticles is considered in sec- 
tion IVllBl while the average interaction in zinc and potassium nanoparticles is discussed in 
section IVII CI Generally, we find that the average phonon- mediated interaction in spherical 
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aluminum and potassium nanoparticles is larger, on average , by a few tens of percent than 
the average phonon-mediated interaction strength extrapolated from the bulk value. For 
zinc particles we find that on average there is almost no deviation from the extrapolation 
from the average bulk phonon-mediated interaction strength. The average phonon-mediated 
interaction in deformed nanoparticles is on average reduced by more than 55% compared 
to the average interaction in spherical particles. The deformations are more effective in 
splitting the HOS and reducing the average interaction in smaller particles than in larger 
ones. 

The addition of Coulomb interaction to the phonon-mediated effective interaction re- 
sults in an overall average attractive interaction in almost all spherical aluminum and zinc 
nanoparticles, while for spherical potassium particles it results in an overall average re- 
pulsive interaction. Most aluminum and zinc particles exhibit an overall average repulsive 
interaction when deformations are taken into account together with Coulomb interaction. 

The energy gap in the three types of nanoparticles is discussed in section IVIIDI Spherical 
aluminum particles exhibit an average energy gap of about 0.13eV when Coulomb interaction 
is ignored. When deformations and Coulomb interaction are taken into account, the average 
energy gap (for all particles for which we were able to find a solution to the gap equations 
with negative condensation energy) is reduced to about 0.025eV. The effect of deformations 
is most pronounced near half-shell filling. The average energy gap in the spherical zinc 
particles without Coulomb interaction is about O.OSeV, while the inclusion of both Coulomb 
repulsion and deformations results in an average gap of 0.007eV. The gap completely vanishes 
in potassium particles when deformations and Coulomb interaction are taken into account. 

Although the average phonon-mediated interaction, when deformations are taken into 
account, is weaker in our system than the bulk extrapolation values, the resulting energy 
gap is still larger than the one found in bulk material because of the single-electron energy 
shell structure. When considering aluminum particles (with both deformations and Coulomb 
interaction taken into account) the average energy gap is larger by two orders of magnitude 
than the energy gap in bulk aluminum. Also, the energy gap obtained in aluminum particles 
(when deformations and Coulomb interaction are taken into account) is typically four times 
larger than the lowest single-electron excitation energy (LSEEE) of the aluminum particles 
and an order of magnitude smaller than the average HOS-LUS spacing. Therefore, the 
energy gap lies in an intermediate energy scale between the LSEEE and the inter-shell 
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Aluminum 


Zinc 


Potassium 


A (GPa)94 


68 


46 


2.8 


^ (GPa)2i 


29 


51 


1.3 


p{gi ■ cm~^)^'* 


2.7 


7.3 


0.9 


Z 


3 


2 


1 


ArF(nm)^^ 


0.049 


0.051 


0.075 


(eV) 


0.033 


0.028 


0.008 


AE^LUS-HOS (eV) 


0.26 


0.33 


0.11 



TABLE I: Basic properties of the three types of nanoparticles con- 
sidered in this work. A and fi are the bulk material Lame constants, 
p is the bulk density, Z is the number of valence electrons per atom, 
Xtf is the Thomas- Fermi screening length, Hlod the average De- 
bye energy of the nanoparticles, and AE'lus-hos is the mean energy 
difference between the HOS and the LUS. The elastic constants and 
density of the aluminum, potassium, and zinc particles correspond to 
T = 0, T = AK, and T = 4:.2K, respectively 

energy spacing. 

B. Average electron-electron interaction in aluminum particles 

The relevant material properties of aluminum, zinc and potassium nanoparticles are sum- 
marized in table [H In Fig. [T] we plot the absolute value of the average effective interaction 
coupling constant G for aluminum nanoparticles, while ignoring Coulomb interaction. The 
coupling constants are plotted for spherical particles and for deformed particles. The in- 
clusion of deformations leads to an average reduction in the average coupling constants of 
about 56%. 

The shell structure is clearly seen in Fig. [H The coupling constant changes abruptly 
when moving from one shell to the next. Small jumps occur in the coupling constant within 
a certain HOS, such as when the number of free electrons in the aluminum nanoparticles is 
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FIG. 1: Absolute value of the averaged phonon- mediated interaction constants of aluminum 
nanoparticles plotted as a function of the number of free electrons in the nanoparticle. The results 
for spherical and deformed particles are shown. The effect of Coulomb interaction is ignored. The 
dashed curve shows an extrapolation from the average phonon-mediated interaction obtained from 
the experimentally measured pairing interaction strength in the bulk, and taking the dimensionless 
average renormalized Coulomb interaction fi* to be equal to 0.1.— 

varied from = 366 to Nf, = 367 or from = 373 to Nf, = 374. These jumps arise from 
an increase in the number of phonon modes in the specific phonon branches that mediate the 
interaction between the electrons in the HOS, due to the increased number of atoms in the 
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FIG. 2: Absolute value of the average phonon-mediated interaction constants —G of aluminum 
nanoparticles containing between 339 to 380 free electrons. All particles belong to the filling of 
the HOS characterized by Z = 10 and n = 1. The jumps in —G between N = 366 and N = 367, 
and between N = 372 and N = 373, are due to the addition of a single phonon to the effective 
interaction. 

particle. Although these jumps may exist, their location and size should not be inferred from 
our calculation, due to the approximate nature of our treatment of the phonon spectrum 
and wave functions. The variation of the average coupling constant when a shell is filled up, 
as well as the small jumps due to the additional interacting phonon, are seen more clearly 
in Fig. [2] where we plot — G for the HOS characterized by / = 10 and n = 1. 
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The effect of deformation is smallest near the opening or closing of a shell. It is more or less 
constant at the intermediate range, especially for the smaller particles. In this intermediate 
range, the energy difference between most electronic levels is larger than the energy of the 
most energetic phonon that can mediate the interaction between the electrons. Thus, most 
matrix elements are in fact equal to zero, and the change in the average coupling constant 
between sequential nanoparticles is small. A more pronounced minimum at the HOS half- 
filling is observed for the larger particles for which the magnitude of the deformation is 
smaller. 

In Fig. in] we plot the aspect ratios of the deformed aluminum particles with even numbers 
of electrons in the HOS and with bulk material elastic constants. The largest deformations 
are obtained for the smallest particles considered, with a maximal aspect-ratio of 1.13. The 
average aspect-ratio for particles with unfilled shells is 1.04. The aspect ratio is calculated 
using the following expression that determines the relation between the amplitude of the 
quadrupole deformation a]^2 [Eq- ( |42l) ] and the aspect-ratio^ 

(64) 



max 



c 2 + 2a,^2 



a 2 - a]!^'' 

where c is the length of the major axis and a is the length of the minor axis. 

We compare our results to an extrapolation from the bulk coupling constant. The phonon- 
mediated interaction coupling constant in bulk material (denoted by G^^) is given by^ 

, _ 4 (A* + /i*) Ej. 
- 3iV^ ' ^^^^ 

where Ef is the Fermi energy of the bulk material. A* is the dimensionless average pair- 
ing interaction strength in the bulk material, and /i* is dimensionless average renormalized 
Coulomb interaction (the dimensionless phonon-mediated interaction in the bulk is usually 
denoted by A). We use the value of A* that is extracted from measurements of the crit- 
ical temperature and the energy gap in bulk aluminum, and for comparison the value of 
A* calculated by Morel and Anderson^^. In order to obtain we use the value /i* = 0.1 
calculated by Morel and Anderson^S for aluminum. The calculation of the phonon-mediated 
interaction by Morel and Anderson is similar in some basic assumptions to ours since they 
used a free electron model, a screened Thomas-Fermi interaction between the electrons and 
the phonons, and they took into account only the effect of longitudinal phonons. How- 
ever, unlike our model, the phonon spectrum in the Morel and Anderson calculation was 
approximated by an Einstein phonon model smeared into a Lorentzian line shape. 
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FIG. 3: Aspect ratio of deformed aluminum particles obtained using bulk aluminum elastic con- 
stants. 

We fit our results to a function of the form 

and list the values of A and a in table HIl The number 750 is just a typical value for the 
number of electrons considered. The values of A for Gl^ [calculated using Eq. ( 165|) ] are 
also given in table [III As can be seen from Eq. (l65l) . a is identically equal to 1 for bulk 
material. We find that G in the spherical nanoparticles is on average larger by about 40% 
than Gl^ extrapolated from the experimental value in the bulk, and by about 15% than 
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Gl extrapolated from the calculations of Morel and Anderson.— When deformations are 
considered G is smaller than G^'^ by about 35% on average. 

Although on average there is no large difference between the smooth extrapolation from 
the bulk interaction and our results, the electronic shell structure causes a large variation 
in the values of G compared to the monotonic behavior of the extrapolation (Fig. [1]). The 
differences between the values of G in the various HOS reflects the differences in the electron- 
phonon matrix elements and the change in the number of phonons that can mediate the 
interaction in a specific HOS. The variation is especially evident when going from a HOS 
with high value of / and low value of ri to a HOS with a low value of / and high value of n 
(or vise versa). For shells with small / we find especially high values of G However, the low 
degeneracy of these shells tends to cancel out the high value of G when the energy gap is 
calculated [Eq. (I58|) ]. 

There is also a large scatter of the detailed matrix elements ^"^^ ^ specific particle 

around the average G of the particle. The large scatter reflects the variation in the electron- 
phonon interaction when considering different m's within a given shell, and especially the 
different number of phonons that contribute to different G^^' 's. The variation in the detailed 
matrix elements increases when deformations are taken into account. 

In a recent work by Croitoru et al.^^ the changes in intra-shell pairing interaction in 
spherical nanoparticles were calculated assuming a uniform interaction, modified from its 
bulk average value only because of the modifications in the electronic wave functions [which 
were taken to be the same as the ones given in equation ([21)]. Croitoru et al.— found a 
large enhancement (compared to the bulk extrapolated value) in the intra-shell interaction. 
However, as discussed above, the effective electron-electron interaction within the HOS is 
anisotropic even within a degenerate HOS, where a large difference between m and m 
results in a reduced G^^i compared to the diagonal matrix elements. The net effect is 
only a moderate increase in the average effective interaction, especially when the HOS is 
characterized by a large value of /. Therefore, we do not expect to find a major enhancement 
in the intra-shell phonon-mediated interaction, especially in large nanoparticles as the ones 
considered by Croitoru et al.— 

Taking Coulomb repulsion into account leads to an overall average attractive interaction 
in the spherical particles as shown in Fig. |H On the other hand, the combination of defor- 
mations and Coulomb interaction results in a repulsive average interaction in about 85% of 
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a 


A 


Spherical - 100% bulk constants 


0.719 


0.0075 


Deformation - 100% bulk constants 


0.757 


0.0030 


Bulk - experimental 


1.0 


0.0057 


Bulk - Morel and Anderson^^ 


1.0 


0.0068 



TABLE II: Fitting parameters for the average phonon-mediated in- 
teraction coupling constants, Gfu = A (^) of spherical and de- 
formed aluminum particles with bulk elastic constants, together with 
the corresponding approximate parameters for bulk aluminum. 



the particles. Both phonon mediated interaction and Coulomb interaction decay with the 
increase in the number of electrons in the particles, but the Coulomb interaction diminishes 
faster. Also, the size of the deviation from spherical symmetry decreases with the increase in 
the size of the nanoparticles. Therefore, we find more nanoparticles with attractive average 
coupling as the size of the particles increases. 

We also evaluate the average effective interaction and energy gap using Lame coeffi- 
cients which are 25% smaller than the corresponding coefficients of bulk aluminum. Using 
these coefficients, the ratio between the longitudinal and transverse speeds of sound remain 
unchanged, and all dimensionless properties of the vibrational modes remain unaltered. 
Therefore, the electron-phonon matrix elements are changed only because of the factor 
cfj = a/p/ (A + 2/i) appearing in Eq. (122|) . The effective electron-electron interaction is 
proportional to the square of the electron-phonon matrix elements, and therefore it is just 
increased by a factor of 4/3 in the spherical particles. 

On the other hand, by decreasing the elastic constants we make the nanoparticles more 
susceptible to deformations. The amplitude of the deformation and therefore the size of the 
energy splitting, is increased by a factor of 1^64727 due to the smaller Lame coefficients [Eq. 
fj42|) ]. Also, the phonon frequencies decrease compared to the particles with bulk elastic 
constants. The lower frequencies together with the larger energy splitting decrease the 
number of phonons that can mediate the interaction between specific electron pairs. The 
net effect is nevertheless an increase of the overall electron-electron interaction, as can be 
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FIG. 4: Average electron-electron interaction of aluminum nanoparticles as a function of the num- 
ber of free electrons in the nanoparticle, taking into account the effect of the Coulomb interaction. 
Results are shown for spherical and deformed particles. 

seen from Fig. [5] and from the fact that a larger portion of the particles examined (31% 
compared to 15%) exhibit an overall average attractive interaction when both deformations 
and Coulomb interaction are taken into account. We note that the ratio between the value 
of G at the closing of a HOS and its value at half filling is larger for particles with reduced 
Lame coefficients than it is for particles with bulk aluminum constants. 
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FIG. 5: Average coupling constants of deformed aluminum particles with Lame coefficients at 75% 
of their bulk value, with Coulomb interaction taken into account. 



C. Average electron-electron interaction of zinc and potassium particles 

The average coupling constants of potassium particles containing 100 to 500 atoms are 
shown in Fig. |6l Since potassium is a non-superconducting metal we can compare our results 
only to theoretical calculations of the strength of the phonon-mediated interaction. The 
calculations of Morel and Anderson yield a major overestimate of the overall dimensionless 
interaction strength. Therefore, we use the results of more detailed calculations^Ii^S, which 
yield values of A (the dimensionless average electron-phonon interaction strength) between 
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FIG. 6: Average coupling constant of spherical potassium nanoparticles without Coulomb repul- 
sion. The extrapolation from the theoretical value of G^^ which corresponds to A = 0.11 is also 
plotted. 

0.11 and 0.16, with most results tending toward the lower values. In Fig. |6] we plot the 
values of G for potassium particles together with the calculated average phonon mediated 
interaction G^^ given in Eq. (165 p . using A = 0.11. The average of the ratio G/G^^ is 
about 1.4 for A = 0.11 and 1.05 for A = 0.16. Unlike the aluminum particles, the addition 
of Coulomb interaction together with deformations results in purely repulsive interaction 
matrix elements. Even when deformation is disregarded, the addition of Coulomb interaction 
results in an overall average repulsive interaction for all potassium nanoparticles. 
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FIG. 7: Average coupling constant of spherical zinc nanoparticles without Coulomb repulsion. We 
also plot the extrapolation from the experimental bulk value, assuming fi* = 0.09.— 

Unlike the potassium particles, the average phonon-mediated interaction in the zinc spher- 
ical particles is larger than the average Coulomb interaction. However, when both defor- 
mations and Coulomb interaction are taken into account we obtain repulsive interaction for 
almost all (99%) zinc particles. The phonon-mediated average interaction of the spherical 
zinc particles is plotted in Fig. [7] and the overall average interaction of the deformed zinc 
particles together with Coulomb interaction is shown in Fig. |8l 
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FIG. 8: Average coupling constant of deformed zinc nanoparticles with Coulomb interaction taken 
into account. 

D. The energy gap 



1. Energy gap of spherical particles 

In Fig. |9] we plot the energy gap for spherical aluminum, zinc, and potassium nanopar- 
ticles, as calculated by the seniority model, without Coulomb interaction. The energy gaps 
are clustered according to the filling of the HOS of the nanoparticles. The gap shows an 
even-odd effect as a function of the number of electrons Nhos in the HOS: A = G{21 + 1) 
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or A = 2GI for even or odd Nhos, respectively. As mentioned in section IVIBt the gap ob- 
tained by the BCS approximation coincides with the seniority gap for even Nhos, regardless 
whether an average interaction or detailed matrix elements are used in the BCS calculation. 
The jumps in G, due to the addition of a single phonon to the effective interaction between 
the electrons, are magnified, by the multiplication of G with half the HOS degeneracy. See, 
for example, in Fig. [9] there is a jump in the gap between Ng = 366 and Ne = 367 or 
Ne = 373 and A^^e = 374 for the aluminum particles; and between Ne = 450 and A^^e = 451 
or A'e = 576 and N^ = 577 for the zinc particles. Figure [TO] provides a magnified view of the 
region A^ = 400 — 500 for aluminum, in order to demonstrate more clearly the jumps in the 
energy gap that are due to transitions into the next electronic shell and that are caused by 
the addition of a phonon to the effective interaction within a given shell. 

The average energy gap of spherical aluminum particles is given in table lllll In calculating 
the average gap for spherical particles we disregard particles with Nhos = 1,4/ + 1,4/ + 2. 
These particles may form a paired state only through the much weaker inter-shell interaction 
which we ignore. Therefore, they possess a smaller or even zero energy gap. 

2. Energy gap of deformed aluminum particles without Coulomb interaction 

The BCS model is used in order to estimate the energy gap for nanoparticles with non- 
degenerate HOS. Therefore, we calculate the energy gap only for particles with even Nhos 
and non-full HOS. The resulting average energy gaps of the aluminum particles for the 
various scenarios considered in this work are summarized in table IIIII In calculating the 
average energy gap of the deformed particles we take into account only particles with even 
Nhos and non-full HOS for which we were able to find a solution to the gap equations 
with negative condensation energy. The fraction of deformed particles (with and without 
Coulomb interaction) with even Nhos and non-full HOS for which we were able to find such 
a solution is shown in table IIIII 

The energy gap is seven times larger on average than the LSEEE of the deformed alu- 
minum particles without Coulomb interaction. In calculating this average ratio (as well 
as the same ratio when Coulomb interaction is taken into account) we take into account 
only particles with even Nhos and non-full shells for which we find a negative condensation 
energy. Of those, we disregard the particles for which the LSEEE involves a shell transition. 
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FIG. 9: (color online) Energy gap of spherical aluminum, zinc, and potassium nanoparticles plotted 

as a function of the number of free electrons in the nanoparticle, without Coulomb interaction, as 

obtained by either the seniority model or BCS model. The difference (equal to G) between the 

gap in particles with an even number of electrons in the HOS and particles with an odd number of 

electrons in the HOS is clearly seen for aluminum and zinc but is hard to observe for the potassium 

nanoparticles due to the relatively small average coupling constant. The larger jumps in the gap 

correspond to a transition from one shell to the next one. The transition between shells is marked 

by three particles with a zero gap, which correspond to Nhos = 1,4/ + 1,4/ + 2. The smaller 

jumps, such as the ones observed for the aluminum particles between A'^e = 366 and = 367 or 

A'^e = 373 and A'^g = 374, or for the zinc particles between Ng = 450 and A'^e = 451 or A'e = 576 

and A'e = 577, correspond to the addition of a phonon to the effective interaction within a given 
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shell. Fig. [10] shows a magnification of the region A'^ = 400 — 500 for aluminum. 
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FIG. 10: Seniority model energy gap of spherical aluminum nanoparticles plotted as a function 

of the number of free electrons in the nanoparticle, without Coulomb interaction, as obtained by 

either the seniority model or BCS model. Five shells are shown - {I = 7,n = 1), {I = 2,n = 4), {I = 

0,n = 5), (/ = ll,n = 1), and {I = 5,n = 3). The transition between shells is marked by three 

consecutive N values for which the energy gap vanishes (within the framework of our model in 

which inter-shell interaction is ignored): Nhos = 4Z + 1, 4Z + 2 for an almost-full or full shell and 

Nhos = 1 for the next shell. Around Ne = 440 one notes 5 consecutive Ne values with a zero 

energy gap. Three of these correspond to Nhos = 4x2 + 1,4x2 + 2, 1 while the other two 

correspond to the shell / = 0, n = 5. The smaller jumps observed between A'^e = 420 and Nf, = 421 

and between = 473 and = 474 are due to the addition of a single phonon to the effective 

interaction between the electrons and are not related to electronic shell effects. 
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Average gap [eV] 


Fraction 


/ gap \ 

\ LSEEE/ 


Spherical - 100% bulk constants 


0.139 


100% 




Spherical - 75% bulk constants 


0.214 


100% 




Deformation - 100% bulk constants 


0.040 


62% 


7.0 


Deformation+Coulomb - 100% bulk constants 


0.024 


51% 


4.1 


Deformation+ Coulomb - 75% bulk constants 


0.045 


50% 


6.6 



TABLE III: Summary of the results we obtained for the energy gap 
of the aluminum particles in the various scenarios considered in this 
work. In the second column, we show the average energy gap of the 
aluminum particles. In the third column, we present the fraction of 
the aluminum particles for which we were able to find a solution to 
the gap equation with negative condensation energy. In the fourth 
column, the average ratio between the energy gap of the deformed 
aluminum particles and the lowest single electron excitation energy 
(LSEEE) is shown. 

In Figs. [TT] and we plot the resulting gap for deformed aluminum nanoparticles without 
Coulomb interaction. We plot the energy gap that corresponds to a solution of the gap 
equations with the lowest negative condensation energy. We note that, for a large portion of 
the deformed particles, a part of the gap parameters Am is essentially equal to zero (within 
the accuracy of our numerical solution). However, this is an artifact that disappears with 
the addition of Coulomb interaction, which results in gap parameters that are both positive 
and negative. 

As can be seen from table IIIII we find solutions to the gap equations with negative con- 
densation energies for most nanoparticles with non-degenerate shells. Also, as can be seen 
from Fig. [11] and more clearly from Fig. [121 the energy gap generally exhibits a quite smooth 
behavior as a given energy shell is been filled up, and a minimum near half-shell fillings. This 
behavior is consistent with the behavior of the average coupling constants (Fig.[T]). However, 
we fail to obtain solutions to the gap equations with negative condensation energy for some 
particles, which many times are located near shell opening or shell closing. These failures 
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are probably due to the limitations of our numerical procedure and not due to a real feature 
of our system. We claim that these are probably failures of our numerical procedure and not 
a real feature of our system because for most energy shells we are able to find appropriate 
solutions for particles located near half-shell filling, where the shell splitting is maximal and 
the formation of a paired state is least likely. We note that although the general structure 
shown in Figs. [11] and [12] seems to be reasonable, we cannot prove that solutions with lower 
condensation energies do not exist. 

3. Energy gap of deformed aluminum particles with Coulomb interaction 

In Figs. [13] and [TH we show the energy gap for aluminum particles when Coulomb in- 
teraction is added on top of deviations from spherical symmetry. The average value of the 
energy gap is given in table lllli We find a solution to the gap equations with negative con- 
densation energy for 50% of the aluminum particles with even Nhos- For these particles, 
we find that the energy gap is approximately 4 times larger on average than the LSEEE. 
However, it should be noted that near shell opening or closing, where the effect of deforma- 
tions is minimal, the gap can be larger than the LSEEE by a factor of 10 or more. Also, 
the average energy gap is still two orders of magnitude larger than the energy gap found 
in superconducting bulk aluminum'^ (which is about O.lSmeV), but smaller (by an order of 
magnitude) than the mean level spacing between the HOS and the LUS (see table [T] for the 
average level spacing between the HOS and the LUS in the aluminum particles ). 

Aluminum particles with smaller Lame coefficients exhibit a larger energy gap on average 
than the aluminum particles with bulk Lame coefficients (see table [In|) . The percentage 
of particles for which we are able to find solutions to the gap equations with negative 
condensation energy is similar to the one obtained for particles with bulk Lame coefficients. 

4-. Gap parameters anisotropy in aluminum particles 

The calculated values of the gap parameters exhibit a large anisotropy compared to 
bulk aluminum. The average ratio between the maximal and minimal absolute values of 
Am is about 10 when both deformations and Coulomb interaction are taken into account. 
By comparison, according to the calculations of Leavens and Carbotte^^, the corresponding 
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FIG. 11: BCS model gap of deformed aluminum nanoparticles, without Coulomb interaction, with 
even Nhos and negative condensation energy plotted as a function of the number of free electrons 
in the particle. The gap is calculated using the interaction matrix elements G^^' , and by solving 
Eqs. (15^1) and (lUTD. 

ratio in bulk superconducting aluminum is about 1.4. We note that Croitoru et al.^^ have 
found a large spatial anisotropy in the gap parameter of spherical nanoparticles, calculated 
by solving the Bogoliubov-de Gennes equations. 
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FIG. 12: BCS model gap of deformed aluminum nanoparticles, without Coulomb interaction, with 
even Nhos varying between Nhos = 562 to Nhos = 612 and negative condensation energy. All 
particles belong to the HOS I = 12, n = 1. The gap is calculated using the interaction matrix 
elements G i, and by solving Eqs. and ([UT|) . 



5. Energy gap of zinc and potassium particles 

In Fig. [15] we show the energy gap of zinc particles when deformations are taken into 
account together with Coulomb repulsion. Although zinc particles are less susceptible to 
deformations than aluminum particles, the smaller electron-phonon interaction and the re- 
sulting effective electron-electron interaction lead to a lower average energy gap for zinc 
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FIG. 13: BCS model gap of aluminum nanoparticles when both deformations and Coulomb repul- 
sion are taken into account plotted as a function of the number of free electrons in the particle. 
Only particles with even Nhos ^-nd negative condensation energy are shown. The gap is calculated 
using the detailed interaction matrix elements G ', and by solving Eqs. and ([UT|) . 



particles than for aluminum particles. We find that the average energy gap in spherical 
zinc particles without Coulomb interaction is equal to O.OSeV, and to 0.007eV when both 
deformations and Coulomb interaction are taken into account. We find a solution to the 
gap equation with negative condensation energy only in 37% of the zinc particles (with even 
Nhos and non-full HOS). We are unable to find solutions to the gap equations for potas- 
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FIG. 14: BCS model gap of aluminum nanoparticles when both deformations and Coulomb repul- 
sion are taken into account plotted as a function of the number of free electrons in the particle. 
Only particles with negative condensation and even Nhos varying between Nhos = 1012 to 
^HOS = 1074 are shown. All particles belong to the HOS / = 15,n = 1. The gap is calculated 
using the detailed interaction matrix elements G^^', and by solving Eqs. (f59l) and (|6T]) . 



sium nanoparticles when both deformations and Coulomb interaction are taken into account. 
This implies that within the framework of the BCS grand-canonical approximation, pairing 
correlations are completely destroyed in potassium nanoparticles. 
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FIG. 15: BCS model gap of zinc nanoparticles, when both deformations and Coulomb repulsion 
are taken into account, plotted as a function of the number of free electrons in the particle. Only 
particles with even Nhos and negative condensation energy are shown. The gap is calculated using 
the detailed interaction matrix elements G^^i , and by solving Eqs. (f59|) and ([6T]) . 

VIII. POSSIBLE MEASURABLE QUANTITIES 

In bulk superconductors, the energy gap is directly observable by measuring the electronic 
DOS of the material (for example, by tunneling experiments), which is significantly altered 
compared to the normal state. By contrast, there is no qualitative difference between the 
spectrum of unpaired electrons in ultra-small nanoparticles and the spectrum of paired 
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electrons, because the spectrum is discrete in both cases. As mentioned above, we anticipate 
that the lowest excitation energy of the paired state [i.e. the energy gap) to be several times 
larger on average than the lowest excitation energy of the unpaired state, and smaller, on 
average, by an order of magnitude than the energy difference between the HOS and the LUS. 
Thus, an observation of excitation energies that lie between the HOS-LUS difference and 
the energy differences expected from the unpaired shell structure may indicate the presence 
of pairing in ultra-small nanoparticles. The effect of the HOS splitting is minimal near 
shell opening or closing, while our results are more accurate near shell closing since the 
number of particles participating in the paired state is maximal. Therefore, like Kresin and 
Ovchinnikov^ we suggest searching for evidence for pairing correlations in nanoparticles 
with Nhos near shell closing. 

Although observation of alteration in the electronic spectrum may be used as an indicator 
of the presence of pairing, one must remember that our description of the single-electronic 
shell structure is quite crude and neglects effects that may contribute to the splitting of the 
HOS, such as non-axially symmetric deformations, surface roughness and disorder. There- 
fore, it is possible that we underestimate the magnitude of the splitting, and as a result 
overestimate the magnitude of the energy gap. If this is indeed the case it will be hard to 
distinguish between splitting caused by pairing and splitting caused by deviations from the 
spherical shell structure. 

The seniority model predicts [Eq. ( !58ll ] that the spacing between adjacent energy levels 
decreases linearly when higher levels within the paired HOS are considered. The presence 
of such a structure in the electronic energy spectrum could indicate the existence of pairing 
interaction. Also, according to the seniority model the difference between two adjacent level 
spacings is equal to 2G {i.e. to a few meV). These two predictions are somewhat modified 
(but not substantially) when the effects of deformations are taken into account within the 
framework of the BCS approximation. 

An alternative fingerprint of the presence of pairing could be found in magnetization 
measurements. Let us assume that a nanoparticle does not exhibit pairing correlations and 
that the degeneracy of the HOS is completely lifted (except for the spin degeneracy) due to 
deformations. This is the accepted description of the energy levels of atomic clusters.— >22 
We further assume that the energy splitting within a sub-shell {i.e. between m and — m) 
is smaller than the splitting between sub-shells {i.e. between different values of \m\), and 
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therefore the electrons first fill up the sub-shells with small values of \m\. We focus on 
nanoparticles with two electrons in the highest populated sub-shell [i.e. half-filling of the sub- 
shell) and consider the zero-temperature value of the particle magnetization. The magnetic 
field inducing the magnetization is taken to be constant, directed in the z direction, and 
weak enough so it does not change the electron population in the various sub-shells. The 
effects of deformations or renormalization on the electronic wave functions are ignored. 

Under the above assumptions we find, to lowest order in the weak magnetic field, that 
the magnetization of the nanoparticle is 

M = -fiBm, (67) 

where /i_B is the Bohr magneton, and we assume that the lowest non-full energy sub-shell 
corresponds to —m. By contrast, if pairing is present in the same nanoparticle, then the 
magnetization will be zero. However, it should be noted that if the number of electrons in 
the particle corresponds to a full highest occupied sub-shell, then the magnetization is zero 
(to the first order in the field strength) also in the absence of pairing. 

A different approach may rely on the detection of the transition between the paired and 
the unpaired state as the temperature is varied. Cao et al.— measured the heat capacity 
of aluminum atomic clusters containing between 43 and 48 atoms. They found a peak in 
the heat capacity of clusters containing 44 and 47 atoms around T = 200K. These peaks 
are claimed to represent the transition between the paired and the unpaired state as the 
temperature of the clusters is varied. The experimental results are somewhat higher, but are 
still in accordance with the calculations of Kresin and Ovchinnikov.— We note that when 
evaluating the energy gap of the nanoparticles at such a temperature, it may be necessary 
to take into account finite temperature broadening of the sharp, almost degenerate, single- 
electron levels. 

IX. SUMMARY 

We investigated the possibility of pairing interaction in metallic nanoparticles containing 
a few hundreds of atoms at zero-temperature. Three materials were considered - aluminum, 
zinc and potassium. We started from a simple model of non-interacting phonons (quantized 
normal modes of vibration of a stress-free elastic sphere) and electrons (fermions in a spher- 
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ical potential box) in the nanoparticles. We introduced an electrostatic screened Thomas- 
Fermi interaction potential between the electrons and the phonons, using the Thomas-Fermi 
screening length in bulk materials. We then derived an effective phonon-mediated electron- 
electron interaction, resulting from the underlying electron-phonon interaction. We found 
that the strongest attractive interaction was between time-reversed electron pairs within the 
same energy shell. We neglected the rest of the effective interaction and assumed a model 
pairing Hamiltonian consisting only of the interaction between electrons in the HOS. 

The effective interaction was derived by means of either the Frohlich transformation^ or 
the similarity renormalization.— The application of the Frohlich transformation was more 
straightforward but restricted to particles with degenerate HOS. However, the degeneracy of 
the HOS is lifted due to deviations from spherical symmetry. The effects of deformations on 
the effective interaction between the electrons were taken into account within the framework 
of the similarity renormalization method. Our application of this technique followed the work 
of Mielke.— However, the interaction we obtained includes a different cutoff function than 
the one obtained by Mielke.— Our cutoff function, which resembles more the result obtained 
by Hubsch and Becker—, reflects the fact that a phonon cannot mediate the interaction 
between electron states with an energy separation larger than its own energy. We also 
included a screened Coulomb interaction in our model Hamiltonian, and found that it plays 
a significant role in reducing pairing effects. 

Unlike the Frohlich interaction, the effective interaction obtained by the application of the 
similarity transformation depends on the renormalized electron energies due to the electron- 
phonon interaction. The size of the renormalization was also obtained using the similarity 
renormalization method. However, within the framework our model, in which inter-shell 
interactions are neglected, the renormalization is unimportant since (for spherical particles) 
it shifts all of the states within the HOS by the same amount. However, renormalization may 
play a more important role when inter-shell interactions are taken into account, especially 
when the non-renormalized HOS and LUS are nearly degenerate. 

On average, there was no large difference between the average phonon-mediated interac- 
tion we obtained and the extrapolation from bulk interaction. However, finite size effects 
induced variation in the values of the average interaction in contrast to the smooth mono- 
tonic behavior of the extrapolation. Finite size effects also caused a large scatter between 
the detailed effective interaction matrix elements for a given particle. This scatter was 
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further enhanced in deformed nanoparticles, and it was important to take it into account 
when we evaluated the energy gap in those particles. We also found that the entrance of an 
additional single phonon to the effective interaction between the electrons resulted in small 
jumps in the average effective interaction, and in more pronounced jumps in the energy gap 
of spherical particles. However, it was hard to observe the fingerprints of this effect in the 
calculated energy gap of deformed nanoparticles. 

The effect of the pairing Hamiltonian on the electronic spectrum was evaluated by using 
the seniority model (when we considered particles with degenerate HOS) or by solving 
the BCS gap equation for the HOS. When Coulomb interaction and deformations were 
ignored, we found that aluminum particles exhibit an average energy gap {i.e. the lowest 
excitation energy of the paired electrons) of about 0.14cV, while the average gap in zinc 
and potassium particles was found to be 0.07eV and 0.02eV respectively. The addition of 
Coulomb interaction, together with deformations, reduced the energy gap by a factor of 5.5 
on average in the aluminum particles and by a factor of 10 in the zinc particles. On the other 
hand, the same effects reduced the energy gap to zero in the potassium particles. Despite 
the large reduction, our calculations indicate that a large portion of both aluminum and zinc 
particles should exhibit pair correlations and modifications in their electronic spectrum. 

We found that within the framework of our approximate model, the resulting energy gap 
is on average intermediate between the LSEEE and the energy difference between the HOS 
and the LUS. Additionally, if pairing is present, then we expect that the magnetization of 
nanoparticles with certain fillings of the HOS (namely, two electrons at the highest sub- 
shell within the HOS) to be different from the magnetization of the unpaired ground-state. 
Therefore, magnetization may serve as an additional fingerprint of the existence of pairing 
correlations in metallic nanoparticles. 

Our results indicate that the size (and maybe even the existence) of the modifications in 
the electronic energy spectrum are sensitive to several factors such as the effective elastic 
constants of the particles, the details of the energy splitting of the HOS (which may be 
affected by additional mechanisms such as disorder, non-axially symmetric deformations, 
and surface roughness), and maybe the detailed description of the Coulomb interaction. 
Therefore, it is probably necessary to consider these aspects more accurately in order to 
reliably predict whether pairing correlations do exist in specific metallic nanoparticles. 

An additional improvement involves applying a canonical-ensemble treatment of the pair- 
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ing Hamiltonian instead of the grand-canonical BCS approximation. The simplest fixed-N 
treatment for the energy gap^*^ is obtained by projecting the BCS ground state, as well as the 
BCS ground state with two "blocked" single electron states, onto the subspace of states with 
a fixed number of electrons in the HOS and calculating the energy difference between them. 
This method yields results that are quite close to the ones obtained by an exact solution of 
the pairing Hamiltonian for aluminum particles with constant interaction strength^iii and 
uniform or random level structure. Therefore, it is probably sufficient in order to evaluate 
the accuracy of the BCS approximation results for our case, where an exact approach is 
inapplicable. More sophisticated methods could also be used, but their application should 
be considered in light of the approximations involved in deriving the effective interaction 
between the electrons and in estimating the single-electron level structure. 

Finally, we note that the standard approach for calculating properties of superconducting 
bulk material relies either on the application of Eliashberg theory^ii^i"— or density func- 
tional theoryi^i^iS^ instead of a BCS approach based on an effective model Hamiltonian. 
In fact, density functional theory was already used to describe pairing in the nucleu o^°^'^^° 
and in nanometric superconductors.- — Specifically, the application of density functional ap- 
proach to equally spaced^i^^ or randomly distributed^ electronic spectrum (both with 
constant pairing interaction) reproduced remarkably well the results of the exact Richardson 
solution. It would be interesting to see how our simple BCS approach would fair compared 
to a density functional approach applied to non-constant pairing interaction and an approx- 
imate electronic shell structure. On the other hand, the adaption of the Eliashberg theory 
to the problem of pairing in ultra-small nanoparticles will enable an alternative (and indeed 
a far more widely used one in the context of bulk material) to our similarity renormalization 
+ BCS solution approach, for the calculation of pairing correlations from the underlying 
electron-phonon interaction. 

Appendix A: Spheroidal phonons and their longitudinal component 

The details of the solution of the linear elasticity equation of motion can be found in 
many references^, here we shall cite only the relevant results. 

The frequencies of the spheroidal modes are determined by solving 

TnT43 - T41T13 = 0, (Al) 
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for 1^0, where 

Til = - ^ - y) MO + 2pRji+iiO , 
Ti3 = /(/ + i)[(/-ih(^)- , 

^43= (^i'-y-^)mv) + vji+M 

where ji is the spherical Bessel function of the order /. For / = one needs to solve 

Til = 0. (A2) 

In the above equations 1] = pR and ^ = ^R, where Cf is the transverse sound velocity. 

The ratio between the amplitude of the longitudinal component A and the amplitude of 
the transverse component C is determined by the stress-free boundary conditions imposed 
on the surface of the sphere {arr = cror = = 0, where aij are the components of the 
stress-tensor). If / > the ratio is given by — Tn/Tis, while if / = then C = and A is 
determined solely by the wave function normalization condition. 



Appendix B: Derivation of the electron-phonon matrix elements 

We derive the explicit expression for the electron-phonon matrix elements Mlmn (EH). 
Using Eqs. ([2]), and ([6])- (Q, we obtain the following expression for the matrix elements 

g-fcTF|ri-r2| 



Mlmn = J j Bi^nJh{Kn^ri)Yil^^{ni)Bi^nJh{h2 



|ri - ral 



^hnsJ ^^''"^"^ 3i3iki3n3r2)Yi^rn3i^2)rldridilirldr2dn2, (Bl) 

where the indexes LMN denote the same as in (fT5l) . We insert the spherical harmonic 
expansion of the screened electrostatic potential (|T3l) . and the explicit expression for Bi^ ([3]) 
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into (IBip . and obtain 

oo I g 

Mlmn = Y.Y1 - 



1=0 m=-l 



R r rri 

drirlii^ {kiri)ji^ (/csri) (kTFri) / dr2rlji^ (/cgrs) ii^ {kTFr2) 

nR 

The angular integrations in flB2p yield 

lY,,^.A^2)Y^U^.)dn2 = 6^^M (B3) 



^ ^^'tlll^+'l)^^ '' (^1, ^2, /; 0, 0, 0) c ih, k, I; mi, -m^, m) (-1)"^^ . (B4) 
We note that the only non- vanishing Clebch-Gordan coefficients in (]B4p are those complying 



with the conservation laws f[T6] - [20|) . Therefore the summation in (IB2p is restricted only to 
these values of /'s, and m's. 

Appendix C: Application of the similarity renormalization 

The similarity renormalization is a renormalization scheme in which an initial non- 
diagonal Hamiltonian, which connects states separated by a large energy difference (up 
to a large initial cutoff A), is transformed, via a continuous set of infinitesimal unitary 
transformations, into a band diagonal effective Hamiltonian with a smaller energy cutoff 
A. The renormalization is carried out perturbatively and in a manner which avoids small- 
energy denominators that may appear in ordinary perturbation expansions, or in a single 
unitary transformation of the initial Hamiltonian such as the Frohlich transformation (for 
more details see^*^). 

Guided by the results obtained using the Frohlich transformation, and in order to simplify 
the notation, we concentrate on the interaction within the HOS and ignore inter-shell inter- 
action. The entire Hamiltonian can be obtained following the same steps we outline below 
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(see also the results of Mielke^^ who derived a general effective electron-electron interaction 
for the bulk system). 

The transformed Hamiltonian H\ can be divided into a diagonal part Hqx and an in- 
teraction part Hi\. It is preferable to carry out the pertubative renormalization expansion 
using only normal ordered operators.— Therefore, we write the Hqx and Hix in the following 
manner 

-^OA = ^ ^hnaX '■ ^l^m^nj^hmzm '■ + ^mA : clngCma '■ (CI) 
iamana ma 

i^IX = ^ ] ^ ] ^ ^ ^mxm2hn3\ ■ '^m\a'^rn2ab\^mi-m2n3 ' ~^ ^ mT.m2hn3^ ' '^\n2a'^'m-\o-^hm\-m2n3 ■ 
mim2 ^3713 a 

+0 (M^) , (C2) 

where nmia is the average occupancy of an electronic state at the HOS, :: denotes normal 
ordering, and the following shortened notation was used 

^m\ ^ImnX 
C-mcr ^Imna 

Mm^rn2hn3X = ^lll^mim2mi-m2X- 

The average occupancy is equal to Ne/{41 + 2) if the HOS is degenerate, where iVe is the 
total number of electrons in the HOS, while in the deformed particles is determined by 
the Fermi-Dirac distribution. The term O (M^) contains new interactions generated by the 
transformation, including the effective electron-electron interaction. 

The renormalization of the matrix elements of a general Hamiltonian matrix element Hijx 
is determined through 

= Uijx [r?A, Hix]i^ + njx^^^^^H.jx, (C3) 
and the matrix elements of the generator of the transformation (denoted by f]ijx) is given by 



where the states i and j are eigenstates of the diagonal part of the renormalized Hamiltonian 
Hqx, with eigenenergies Eix and Ejx- The properties of the function Uijx are discussed below. 
The detailed derivation of Eqs. (]C3|) and (]C4p can be found, for example, in the original 
paper by Glazek and Wilsoii^ or in the work of Mielke.— 
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The function IlijX IS cL continuous and differentiable function of the argument 
(A — \Eix — Ejx\) X (3. It is equal to one if (A — \Eix — Ejx\) x f3 « and to zero if 
(A — \Eix — Ejx\) X 13 » 0, where (3~^ is the width in which Uijx varies from one to zero 
around \Eix — Ejx\ = A. We use a single value of (3 to all Uijx and take (3~^ to be small so Uijx 
varies sharply around \Eix — Ejx\ = A. The meaning of "small" and "sharp variation" 
is given in what follows. 

The function Vijx is defined as Vijx = 1 — Uijx- The factor Vijx multiplying the right-hand- 
side (RHS) of Eq. f lC4|) [and thus also the RHS of Eq. ( 1C3|) ] ensures that energy denominators 
small compared to A {Eix — Ejx < A) are avoided. 

Eqs. (1C3P and (IC4p cannot be solved explicitly. Instead, we assume that it is possible 
to expand the generator rjx and the transformed Hamiltonian Hx in powers of the electron- 
phonon coupling coefficients Mmirn2hn-i\ and disregard any terms in Hx which are propor- 
tional to the third power of Mmirn2hn3\ or more. 

The lowest order contribution to rjx comes from the second term in the parentheses in Eq. 
(IC4I) and is linear in Mmrm.2hhn3>^- The factor rijx ensures that Hqx do not contribute to rjx- 
The renormalization of Hx is affected by the generator through the commutator on the RHS 
of Eq. flC3|) . Thus, only the lowest order term of the generator affects the renormalization of 
the Hamiltonian up to the second order in Mmim2i3n3X- We obtain the following expression 
for rjx 

V>. ~ ^ ^ ^ ^ ^ ^ Vmim2l3n3\ ■ (^mia^rn2<jbl3mi-m2n3 ■ ~Vmim2l3n3\ ■ ^m2cr^'m-\crbl3n3m\-m2n3 ■ 
mim2 l3n3 cr 

+0 (M^) , (C5) 



where 



and 



_ rmim2hn3\ dhlUmim2hn3\ \ D ( A/T^\ ( 

Vmim2l3n3X — TT ^Wmimaisns + <->' j , (^Ubj 



Mmim2«3n3A — — j^miA — £m2X + f^hn^x]) X /?] • (C7) 

We note that we do not take into account corrections to the electron-phonon matrix 
elements arising from deviations from spherical symmetry. Therefore, the electron-phonon 
interaction cannot change the energy structure within the HOS (although the inter-shell 
energy structure is modified) and all electronic energy differences e^iA — £m2X are in fact 
A- independent. Thus, in the rest of the calculation, as well as in the expression for the 
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renormalization of the generator matrix elements (lC6p . we replace EmiX — £m2X with the 
initial energy differences £miA — £m2A, where EmA are the non-renormalized energies of the 
electrons. 

Inserting the generator fICSp into Eq. flC3p we find that the commutator in Eq. fICSp 
contributes to the renormalization of Mmim2hn3X only to the third power (or more) of the 
electron-phonon coupling. Therefore, the renormalization is determined through the follow- 
ing equation 

-J^ — rmim2hn3X jWmimaiansA + '-^ [^^^ ) i l*^oj 

with the solution 

■^mim2l3n3X ■^m\m2l3n3K^mim2l3n3X -^m\m2l3n3\'^m\m2l3n3X(i 1 2 3 3. (^9) 

Note that the function emim2«3n3A has the same asymptotic behavior as Um^m2i3n3X- 

The renormalization of the single-particle electron and phonon energies is determined 
solely by the commutator in Eq. (IC3p . while the second term in fICSP is irrelevant due to the 
presence of the pre-factor rm^rn2i3n3X- By contrast, both terms contribute to the generated 
electron-electron interaction. We calculate the commutator and use the expressions for 
the matrix elements of the generator flC6P and the renormalized electron-phonon coupling 
coefficients (]C9p . in order to derive differential equations describing the renormalization of 
the single-particle energies Smux and w/gnaA, and the effective interaction V^/jv/'a between 
time-reversed electron pairs {{M t, -M ;} and {M f, -M ;}). We obtain the following 
three differential equations 



'^^MM'X _ „ ( ^ ^ '^^MM'l3n3X I 

J 1 p '^^M'M;3n3A I I 2 

^ , ^ I fc, , ^MM'l3n3X Tx V^^MM'l3n3A\ 

^M'A " ^MA + iT^hniX "A 

I '^MM'xdUMM'X-tr (n'\r\\ 

+- ;tt — Vmm'x^ (^10) 

"'MM'X 



deMsx _ \- [ 1 - n.m,„ 2 de^Mi: 

— 1\ ~ I l^^-^mMhnsAl TT" 

^ — V^mcrA — £^MEA + i^hniX "A 
mt3n3cr ^ 

de" 



^ma I ,^ |2 "^MmiansA \ ^m^\ 

IMMrnhnsAl TT , {^^^> 

^MSA — £maA + n^lsnSX / 

duJl^nsX \ ^ Tl'inia '^mi-m3(T i ,^ |2 '^^mimi-m3Z3n3A /pi r,\ 

TT ~ / ^ ~ ~T M^^mimi-m3i3n3A| TT k^^^) 

"A , ^rriiuA tmi— mscrA + "^/srtsA wA 
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where 

^MM'X = u[{\-2 |£A/aA - ^mVaI) >< ■ (C13) 



From Eq. (1C12|) we see that the renormahzation of the phonon frequencies is zero (at least 



up to the second order in Mmim2hn-i\) when we consider spherical nanoparticles with degen- 
erate energy shells. Even if deviations from spherical symmetry are taken into account, the 
corrections are small compared to the non-renormalized frequencies. Therefore, we neglect 
the renormahzation of the phonon frequencies and assume that they are A- independent. 
The solution of (ICIOI) is given by 

VWm'a_ = -exp ( / p{t)dt \ I exp (- j p(t)dt\ g{s)ds, (CM) 
where we used the fact that V^/a^'a = 0, and the functions p{X) and g{\) are defined as 

jx\ _ ( \^ 1 ^^MM'ljnsX I |2 

£MaA - ^M'aA + '^hn'iX CJA 

(E^U + E^J, (C15) 

\l3n3 h^i / 

p(A) = _ ''MM'xdUMM'X _ (C^g) 

'^mm'x dX 

In order to explicitly calculate Vj^jj^^'^^,^ we need to choose a specific form of Uj^j^^'i^^^^ 
and Uj^jj^iy^. We define rj^j^'i^^^^ and J^a/a/'a Fermi-Dirac functions with a width that 

we denoted before as f3~^. These functions (as well as Uj^j^^^i^^^^ and Uj^,^^.f^) vary sharply 
from zero to one around A = X^^^'i^^^ and A = X^.jm' respectively, where X^^'i^n^ and A^^^' 
are defined as 

^MM'hm = l^MaA- £M'aA + f^^hnsl, (C17) 



'^mm'x 



Each term in the two sums in flClSP is characterized by certain A^.^j^^';^^^ and A^j'j^^;^^^ at 
which the functions ^mm'^usX ^m' Mhri^x ^^^y from one to zero as A is lowered. These 
variations result in formation of peaks in g{X) around the various AA/Af'igna X^i j^i^^^. 
We take to be much smaller than any energy difference in our system, and therefore 
ensure that the variation from one to zero of ^MM'hnzX ^"^^ ^M'MismX ^^^^ enough, so 
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there is only a small overlap between the various peaks in g (A). The smallness of also 
ensures that there is only a small overlap between the peak located around the lowest value 
of ^MM'hns drop (from one to zero) in p(A) at X^m' ■ 

If effects of deformations on the electron-phonon matrix elements are taken into account, 
then the difference Sm^x — e^-iX becomes A-dependent. In this case we need to impose an 
additional condition on the size of We choose to be small enough so the change in 
\eMa\ - ^M'axl ^s A is Varied across ^MM'hn^ ^mm'^ is small compared to either XMu'hni 



or Xmm'- This last condition can be written as 



d \£M(t\ — £ 



M'aW 



dX 



-1 

(C19) 



for all A^^jv/' and XMM'hn:,- 

Using the specific form we chose for Uj^i ^^ji^^.^)^, we can show that ej^i j^^^^^y^ is essentially 
equal to w^'M/snsA even in the transition zone around XMM'ian^ (and not just as A tends to 
zero or to A). We can also show that the terms in the first sum of ^'(A) fIClSI) with XMM'hna 
larger than both Xm'misus and X^m' a finite non-zero contribution to V^M'A„i„' as long 
as Xmin is smaller than Xf^j^fi^^^ but larger than Xp^^f. Furthermore, due to the narrowness 
of the variation in gl^^ (A) around X^^j^'i^n, the integral ft^^^exp J^^ p{t)dt^ Q^n-Mds 
is independent of A^m, as long as Xmin is slightly smaller than Ajy,^^^';^^^. 

On the other hand, all terms in the first sum of g{X) with X^jj^^i^^^^ smaller than either 
'^m' Mhns -^A/Af' Si'^^ exponentially small contributions, regardless of the exact value of 
Xmin- The same is true for the second sum in ( ICISP with the roles of Aj^j^'^^^g and Xj^ji j^i^^^ 
interchanged. 

As long as Xmin is larger than Xj^^^i the exponential pre-factor in ( ]C14I) is approximately 
equal to one. However, if Xmin is chosen to be smaller than A^^j^^' then the pre-factor becomes 
exp (— A^^j^^'/3), and the effective interaction is suppressed. This is in accordance with the 
general scheme of the similarity renormalization, since Xj^j^f is just the difference between 
the single-electron energies of the pair {M '|", —M J,} and the pair {M' '|", — M' 4}- 

Therefore, if we try to completely diagonalize the Hamiltonian (up to the second order 
in Mm^m2hnz\) by taking Xmin to zero, we obtain decoupled free electrons and phonons with 
renormalized energies. This is not surprising since our application of the similarity renor- 
malization method is perturbative in its nature, and one cannot hope to obtain the paired 
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ground-state of the electrons from the single-electron ground-state via a perturbation expan- 
sion. However, our aim is not to completely diagonalize the Hamiltonian but to decouple 
the electrons from the phonons sufficiently so as to obtain the effective electron-electron 
interaction while accounting for all relevant mediating phonons. This goal is achieved by 
carrying the integration in Eq. (]C14p down to Xmin which is slightly smaller than the smallest 
of >'MM'hn, or Xm'misu, that is still larger than X^j^>. 

The contribution from the first sum in g{X) flC15|) to V"^/m'a„,„ (^^ long as Xmin is larger 
than Ajyj^' and smaller than the smallest relevant X^^j^fi,^,^^ and A^j'^;^^^) is 

O I fij 1 2 

TAfirst sum _ ^ V^-^MM'hn^M (C9C\\ 

hn-i 

where Smk > ^m'a ^ma — ^m'a < ^hn-i- The contribution from the second sum is 

TAsecond sum _ \ ^ ^ \-^MM' l^nsA] (C'9^ \ 

hns M'A ^MA + il^hns 

where Sj^f^ > Sma and ^j^/'y^^ — Sma < ^i-^na- Both contributions do not depend explicitly 
on Xmin- We see that either the first or the second term contribute to the interaction matrix 
element, depending on the sign of Ema — ^m'a- Therefore, we can write the interaction 
matrix element as 



^MM' = AfM fgnaA Q [hWi^^^ - {EmA - £m'a\) ■ (^22) 

^ \£ma - £m'a\ + l^hna 

Inspecting the expression in (]C22|) we can see that the resulting effective interaction is 
always attractive and decays monotonically with the increasing separation between the elec- 
tronic energy levels. Furthermore, a phonon cannot couple electron pairs with energies that 
are separated by more than its energy. The result in ( 10221) coincides with the effective inter- 
action obtained by the application of the Frohlich transformation f l33|) . if the nanoparticle is 
completely spherical. We note that if inter-shell interaction is not neglected then Eq. ( ]C22p 
includes a summation over all relevant shells as well as inter-shell terms. In this case the 
non-renormalized energies of the electrons in Eq. (IC22p are replaced by the A-dependent 
renormalized energies. 

The renormalization of the single-particle electron energies is obtained by integrating Eq. 
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(ICllI) . The integration from Xmin up to A yields 



I I^MMmZansM (C23) 

The sum ^ in ( ]C23P runs over all values of EmaA and huji^n3 for which XmMhm is larger 
than Xmin, and the sum 2J runs over all values of EmaA and hwi^ns ^oi which XMmhns is 
larger than Xmm- In order to take into account the effect of all relevant phonons we need to 
take Xmin to be smaller than the energy of the least energetic spheroidal phonon that can 
interact with electrons in the HOS. 

We now consider an initial untransformed Hamiltonian that includes a Coulomb inter- 
action term. In principal, the additional Coulomb term should modify the generator and 
the flow of the Hamiltonian. However, if we retain the generator (ICSP and do not include 
in it terms arising from Coulomb interaction, then the flow of the electron-electron effec- 
tive interaction remains unmodified up to the second order in the interaction coefficients. 
Therefore, the only modification introduced into the calculation is the appearance of the 
term M^^^, on the RHS of Eq. flC22p . since the electron-electron interaction is equal to the 
Coulomb term and not to zero at A = A. Furthermore, the inclusion of the Coulomb term 
in Hjx does not affect the renormalization of the electron energies as long as the generator 
(1C5I) is retained, since the commutator on RHS of (1C3I) do not yield additional terms that 
contribute to the renormalization of the electron energies up to the second order (including) 
in the electron-phonon and Coulomb interaction coefficients. The generator (1C5I) cannot 
induce the elimination of the Coulomb interaction, which even after the transformation is 
still able to couple states with large energy difference. 

We note that Mielke^ solved the equivalent Eq. to ( IClOP in the bulk while assuming 
'^kk'qx — 1' which is equivalent to assuming w^y^-jv/'A = 1 in our calculation. This assumption 
leads to p(A) = and to the disappearance of the factor m^^^/';^, which multiplies the sum 
in ( IClSp . This last factor is responsible for the Heaviside function in ( ]C22[) . Therefore, in 
Mielke's interaction, phonons with energy smaller than \eMA ~ ^m' a\ "^^^ contribute to the 
matrix element V^^j' and therefore artificially enhance the effective interaction. Further- 
more, the effective interaction of Mielke is not eliminated in the limit A — t- 0. This result is 
inconsistent with the general scheme of the similarity renormalization method. 
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It should be noted that Mielke investigated a model in which the effective interaction was 
mediated by non-dispersive phonons with an energy that is much larger than the energy of 
the electrons. In this model, and as long as Amm is taken to be slightly smaller than the 
single frequency of the non-dispersive phonons, our results and the results of Mielke coincide. 
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